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ABSTRACT 


Maximus  likelihood  methods  are  applied  to  a  series  of  aonopulse 
problems,  involving  both  angle  estimation  and  signal  detection.  Only  the 
two-beam,  off-boresight  monopulse  problem  is  studied.  Explicit  maximum 
likelihood  estimators  are  obtained  in  Part  I,  and  their  probability 
distributions  will  be  discussed  in  the  forthcoming  Part  II.  Both 
deterministic  and  stochastic  signal  models  are  used  here,  and  maximum 
likelihood  estimates  are  obtained  for  the  single  pulse  case  and  for  different 
models  of  the  multiple  pulse  problem.  Particular  emphasis  is  given  to  the 
problem  of  angle  estimation  in  correlated  noise,  representing  the  case  of 
arbitrary  Interference. 
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PART  I:  STRUCTURE  OF  THE  ESTIMATORS 


1.  INTRODUCTION 

In  classical  nonopulse ,  a  receiver  has  two  matched  channels,  each 
connected  to  its  own  antenna.  The  antennas  have  different  patterns,  or 
different  phase  centers,  and  signals  of  known  waveform  arrive  from  an  unknown 
direction  which  is  to  be  measured  by  the  receiver  system.  This  entire  study 
is  directed  at  the  rather  special  case  in  which  the  signal  direction  is 
characterized  by  a  single  unknown  parameter,  such  as  azimuth  or  elevation. 

This  restricted  problem  is  interesting  because  of  its  long  history  and  also 
because  it  lends  Itself  to  a  fairly  complete  mathematical  analysis  in  an 
Interesting  way. 

The  results  developed  here  have  practical  application  to  direction 
finding  systems  in  which  only  one  component  of  direction  is  unknown,  and  also 
to  the  case  where  the  antenna  gains  are  assumed  insensitive  to  the  other 
component,  as  with  a  radar  using  fan  beams.  In  systems  measuring  both  azimuth 
and  elevation,  the  techniques  discussed  here  can  be  used  separately  on  each 
component,  with  somewhat  suboptimum  results. 

The  several  eatlmatlon  problems  analysed  here  are  modeled  and  solved  in 
rather  general  terms.  However,  the  baelc  monopulse  problems  which  motivated 
the  work  are  spelled  out  in  enough  detail  to  keep  the  results  in  direct  touch 
with  practical  applications.  The  method  of  maximum  likelihood  is  used 
throughout,  and  the  focus  is  on  fixed-sample  angle  estimation  (i.e., 
off- bore-sight  monopulse),  and  not  on  angle  tracking,  as  such. 

A  mathematical  setting  has  been  developed  for  the  special  problem  of 
two-beam,  one-parameter  monopulse,  which  appears  to  be  unique  to  this  case. 
Instead  of  dealing  with  bean  gain  ratios  and  complex  data  sample  ratios  as 
points  in  a  complex  plane,  this  plane  is  stereographlcally  mapped  onto  the 
unit  sphere.  The  likelihood  function  is  directly  related  to  distance  on  the 
surface  of  this  sphere,  so  that  maximum  likelihood  estimates  have  very  simple 
interpretations  as  geometrical  projections.  Transformations  of  coordinates, 
of  which  the  well-known  transformation  of  amplitude-comparison  to 
phase-comparison  monopulse  is  an  example,  correspond  to  rigid  rotations  of  the 
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sphere,  and  a  slightly  more  general  class  of  transformations  allows  for  the 
"whitening"  of  problems  Involving  external  interference  of  any  kind. 

The  remaining  sections  of  Part  I  deal  first  with  the  simple  case  of  one 
sample  pair  and  white  noise,  where  the  spherical  model  is  introduced  and 
illustrated  in  terms  of  standard  amplitude-comparison  and  phase-comparison 
systems.  Transformations  are  discussed  next,  and  it  is  shown  that  the 
amplitude-comparison  and  phase-comparison  formulations  are  different  versions 
of  the  same  general  problem,  viewed  in  special  coordinate  systems. 
Generalizations  are  then  introduced  along  two  lines.  In  one  line,  multiple 
sample-pairs  are  considered,  with  coherent  and  incoherent  signal  models.  In 
the  other  line,  non-white  noise  is  introduced,  permitting  the  analysis  of 
direction-finding  in  the  presence  of  Interference.  Finally,  a 
non-deterministlc  signal  model  is  used  to  show  that  the  same  estimators  for 
direction  of  arrival  are  obtained  as  in  the  deterministic  models  of  the 
previous  sections. 

In  Part  II,  to  be  published  separately,  the  performance  of  the  maximum 
likelihood  estimators  is  discussed.  Explicit  probability  density  functions 
are  given  there  for  the  single  pulse  case,  for  both  Internal  and  external 
sources  of  noise.  For  the  multiple  pulse  case,  a  recursive  expression  for  the 
probability  density  function  is  given,  together  with  explicit  formulas  for 
certain  moments  which  can  be  used  to  characterize  accuracy. 

Many  of  the  results  in  this  report  have  been  obtained  by,  or  in 
collaboration  with,  J.  R.  Johnson.  His  insight  and  encouragement  has  played 
an  essential  role  in  the  development  of  the  ideas  presented  here. 
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2.  GENERAL  FORMULATION  OF  THE  CLASSICAL  TWO-BEAM  MONOPULSE  PROBLEM 
In  a  two-channel  monopulse  system  we  model  the  complex  modulation 
functions  at  the  respective  antenna  terminals  as  Zj(t)  and  Z2(t),  where 

Zk(t)  s  Aeifi  S(t)Vk(Y>  +  Wk(t)  , 

for  k  -  1,2.  The  rf  amplitude  and  phase  are  represented  by  A  and  6,  while 
S(t)  stands  for  the  complex  signal  waveform,  assumed  known.  Azimuth  (or 
elevation)  is  called  y,  and  the  Vk(y)  are  the  complex  voltage  gains  of  the  two 
antennas  in  direction  y.  Finally,  W^(t)  and  W^Ct)  are  totally  independent, 
complex  circular  white  noise  processes,  each  with  zero  mean  and  single-sided 
power  spectral  density  N0.  This  last  assumption  results  in  the  complex 
covariance  function 

E  Wk(t)Wk*(f)  -  2NC  «(t’-t) 

for  each  noise  process,  and  the  circular  property  is  represented  by  the 
equation 

E  Wk(t)Wk(t')  -  0 

for  k  ■  1  or  2. 

These  signals  are  passed  through  identical  filters,  matched  to  the 
expected  modulation,  so  that  the  complex  outputs,  sampled  at  the  correct  time 
for  maximum  signal  components,  are 

zk  -  G  /  S*(t)Zk(t)dt  =  sk  +  nk  . 

The  real  and  imaginary  parts  of  each  complex  sample  represent  in-phase  and 
quadrature  components,  and  G  is  a  gain  factor.  The  separate  signal  and  noise 
terms  are,  of  course, 


sk  -  G  Aeifi  Vk<Y>  J|S(t)|2  dt 


nk  -  G  /  S*(t)Wk<t)dt  , 


and  signal  parameters  have  been  suppressed  in  the  notation. 

The  noise  samples  are  complex,  circular  Gaussian  variables,  and  we  have 
the  expression 

Ejnk|2  -  2G2  N0  /  |S(t)|2dt 

for  the  variances  (ksl,2). 

We  now  choose  the  filter  gain  factor  to  make 

g2n0  /|S(t)|2  dt  -  1  , 

so  that 


E I nk |  -  2  . 

Thus  the  real  and  Imaginary  parts  of  n  and  n  constitute  four  real  Gaussian 
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variables,  each  with  zero  mean  and  unit  variance.  Next,  we  introduce  the 
normalized  voltage  gains 


vk(Y)  = 


VT) 


|Vi(T)|2  +  |v2(y)|2 


and  write  the  signal  components  in  the  form 


8k  “  *>vk(y)  , 

where  b  is  our  last  new  parameter,  representing  complex  amplitude,  and  given 


b  -  G  Aei6\/|VI(Y)|  +  |V2(y)|  /|S(t)|  dt  . 
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The  argument  of  b  is  still  the  rf  phase ,  and  the  magnitude  is  determined  from 


M2  -  §£■  i|Vi(Y)|  +  |v2(y)|  }  /|s(t)|  dt  , 


as  a  result  of  our  choice  of  G.  But 


|  A2|Vk(Y)|2  /|s(t)|2dt 


is  the  total  signal  energy  in  the  kth  channel  at  the  antenna  terminals,  and  if 
the  sum  of  these  two  energies  is  called  E8,  the  total  signal  energy  collected 
by  the  system,  then  we  have 


and  the  signal  to  noise  ratio  of  the  system  may  be  defined  as 


SNR  - 


l«ll‘+  l«2^ 

E|ni|  +E|n2| 


Everything  that  follows  is  based  on  the  simple  model 


*k  -  bvk(y)  +  nk 


with  normalised  beam  gains  and  univariate  noise  components,  as  defined  above. 
It  should  be  clear  that  the  special  assumptions  of  ideal  matched  filtering  and 
sampling  can  be  relaxed ,  with  only  a  change  in  the  significance  of  b.  In 
Part  I  we  are  concerned  only  with  the  estimation  of  y ,  and  b  is  a  nuisance 
parameter,  but  the  connection  with  signal  energy  and  noise  spectral  density  is 
important  in  Part  II,  where  performance  is  evaluated  and  signal  to  noise  ratio 
is  a  fundamental  quantity.  It  is  then  important  to  relate  b  to  physical 
parameters  by  means  of  a  specific  model,  as  we  have  just  illustrated  above. 


yr*. 


Note  that  no  aasumptions  have  been  Bade  with  respect  to  £he  voltage  gains uof 
the  antennnas,  and  that  the  basic  model  could  apply  to  any  problem  jjj^ 
characterized  by  sample  pairs  with  relative  signal  components  fixed  by  a 
s Ingle  real  parameter. 

The  probability  density  function  for  the  two  complex  samples  is 


where 


f(*l»*2)  “  exP  )  * 

4w  L 


L2  -  |zx  -  bvx<Y>|2  +  |z2  -  bv2(Y)|2 


Maximum  likelihood  (ML)  estimation  of  b  and  y  is,  of  course,  equivalent  to 
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finding  the  values  which  minimize  L  .  When  we  expand  this  expression,  using 
the  normalization  property  of  the  v^,  and  complete  the  square  we  find 


L2  -  |*i  | 2  +  |z2|2  +  |b  -  tv1*(y)z1  +v2*(y)z2]|2 


-  |v1*(y)«i  +  v2*(y)z2| 


The  estimate  of  b,  given  y,  is  of  course 


b(Y)  =  vx*(y)si  +  v2*(y)*2  * 


and  y  must  minimize 


Min  L  -  |zx|  +  | *2 |  -  |vx*(y)zi  +  v2*(y)z2|  , 

D 


after  which  &  ■  b  (y)  « 

Again  we  expand ,  and  write 


Min  L 
b 


2  2  |v1(y)|2|*i|2+|v2(y)|2|z2|  +  2Re[vi*(Y)v2(Y)2i22*] 

(1*1  +  I*  I  ){1 - - - - - > 

l«ll  +  I z2 I 


The  quantity  in  curly  brackets  will  determine  the  estimator  of  y>  and  it  is 
clear  that  this  expression  depends  on  the  data  samples  only  through  the 
complex  ratio  Z}/Z2  *  Recalling  the  normalization,  it  is  also  obvious  that 
this  quantity  depends  on  the  gain  factors  through  the  ratio  v^/v2  only. 

Ue  express  these  ratios  in  terms  of  four,  real  angle  variables,  as  follows: 

vi(y) 

— — r  5  tan(0/2)e*$ 


—  =  tan(n/2)eM> 

*2 


The  angles  6  and  n  are  restricted  to  the  range  [0,x],  so  that  4  and  4)  are  the 

arguments  of  the  respective  ratios*  Angles  6  and  4  depend  implicitly  on  y>  of 

course,  and  we  may  interpret  8  and  t|  as  polar  angles,  <j>  and  4  as  azimuthal 

angles,  describing  points  on  the  unit  sphere.  The  data  ratio  is  mapped  into  a 

single  point  (n,4)»  but  the  gain  ratio  traces  out  some  trajectory  on  the 

sphere,  as  y  varies  over  its  normal  range. 

The  direction  finding  properties  of  the  actual  system  are  completely 

specified  by  the  "characteristic  trajectory”  over  which  (0,4)  varies,  and 

there  is  no  guarantee,  in  general,  that  this  trajectory  is  a  simple  curve.  It 

is  also  not  guaranteed  that  any  point  on  the  trajectory  corresponds  to  a 

2 

8 ingle  value  of  y ,  but  it  is  clear  from  the  form  of  L  that  there  will  be  no 
way  to  resolve  these  ambiguities  from  the  data  samples,  z1  and  z2> 


7 


It  is  easy  to  show  that 


Iv^y))2  -  sin2(6/2) 
|v2(y)|  ■  cos  (0/2) 


|  (1  -  CO80)  , 

|  (1  +  co80),  and  that 


vi(y)v2*(y)  ■  sin(0/2 )cos(6/2 Je1*  -  r  sln0  e1^ 


Similarly, 


1*1 1 


l«il2  ♦  M2 


_  1 


j  (1  -  cosn)  , 


1*2 1 


1*1 1  +  l*2l 


■  j  (l  +  cosn)  ,  and 


*1*2* 


I *1 1  +  l*2l 


■  r  Sinn  * 


2  2 


When  these  expressions  are  substituted  we  obtain 

Min  L2  ■  l(|*i|2+|*2|2)tl  -  co80  cosn  “  sin0  sinn  cos($  -  ♦)]. 
b  ‘ 


But 


cos6  cosn  +  sin8  sinn  cos($  -  <i)  ■  cos  A  , 


where  A  is  the  arc  distance  between  (0,$)  and  (n,i|>)  on  the  surface  of  the 
sphere.  Finally, 
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.  A'J. 


Min  L2  -  ( |*i t 2  +  |*2 1 2>  sin2 (A/2)  , 

b 

and  the  estimation  problem  reduces  to  minimizing  A,  or  finding  that  point  on 
the  characteristic  trajectory  of  the  system  which  is  closest  (in  the  ordinary 
metric  sense  on  the  sphere)  to  the  data  point  (n,<|))  . 

If  the  characteristic  trajectory  is  complicated ,  then  the  spherical 
interpretation  is  of  no  particular  help  in  finding  an  algorithm  for  explicit 
minimization,  but  in  some  important  idealizations,  it  permits  a  simple, 
explicit  answer  to  the  estimation  problem. 

The  equation 

w  -  tan(6/2)ei* 

establishes  a  mapping  of  the  complex  w-plane  onto  the  unit  sphere  and  vice 
versa.  It  is,  in  fact,  the  well-known  stereographic  projection  of  the  sphere 
onto  the  plane.  The  plane  may  be  taken  as  the  extension  of  the  sphere's 
equatorial  plane,  with  points  geometrically  projected  between  sphere  and 
plane  along  lines  originating  at  the  south  pole.  The  mapping  is  conformal  and 
also  maps  circles  into  circles.  We  do  not  need  any  of  these  facts  (although 
they  help  the  intuition) ,  since  it  is  more  fundamental  to  this  problem  to 
concentrate  on  the  many-to-one  mapping  established  between  complex 
two- component  vectors,  such  as 


and  the  corresponding  points  (n,^)  on  the  sphere. 

In  addition  to  the  estimation  of  signal  parameters,  it  is  often  necessary 
to  decide  on  the  presence  or  absence  of  any  signal  in  the  first  place.  The 
"noise-alone "  hypothesis,  b*0,  is  characterized  by  the  sample  probability 
density  function 
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) 


where 


l  L02 

“  — 9  **P  (”  “5 — 
4w  1 


t 


L02  -  |*i|2  +  |*2|2  • 


Maximum-likelihood  testing  then  reduces  to  evaluation  of  the  test  statistic 

2  2  a  a 

LO  -  girj  L  .  If  the  ML  estimate,  y ,  corresponds  to  the  minimum  distance  A  , 
then 


Min  L2  -  r  < | *i | 2  +  |*i|2)(l  -  cosA  ), 
b,Y 


and  detection  is  based  on  the  test 

j  ( | z 1 | 2  +  | Z£ | 2 ) ( 1  +  cos  A)  >  threshold  . 


Examples  of  the  ML  estimation  of  y  and  the  corresponding  detection  criteria 
are  given  in  the  next  two  sections,  beginning  with  the  familiar  examples  of 
phase-comparison  and  amplitude-comparison  monopulse. 

Before  leaving  the  topic  of  detection,  it  is  interesting  to  consider 
a  refinement  sometimes  called  interference  detection,  or  "data  editing". 
Monopulse  measurements  are  easily  disturbed  by  extraneous  signals,  since  the 
estimation  procedure  depends  critically  on  the  assumption  of  a  single  source, 
with  the  consequent  cancellation  of  the  unknown  signal  amplitude  and  phase  in 
the  ratio,  z^z2  *  uP°n  which  estimation  is  based.  The  presence  of  another 
signal  component  can  be  detected,  however,  and  this  fact  used  as  a  flag  to 
inhibit  processing  of  the  affected  data  for  signal  location. 

Interference  detection  is  accomplished  by  adding  a  third  hypothesis, 
"multiple  targets",  to  the  single-target  and  noise-only  hypotheses  already 
discussed.  If  two  signals  are  present,  the  pdf  of  the  sample  pair  has  the 
standard  form 
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f(*l.  *2> 


9 


_L  -<l2/2) 

4x2 

which  applies  to  all  the  hypotheses,  but  in  this  case  L2  will  be 

2  2  2 
L  ■  |*i  -  bv !  -  b'y/l  +  \z2  ~  bv2  -  b'v2'|  . 


Here  b  and  b'  are  complex  amplitudes,  vl  and  v2  are  voltage  gains  for  the 
first  signal  component,  while  vl'  and  v2'  refer  to  the  other.  The  second  term 
can  be  made  to  vanish  by  the  choice 

b*  -  (*2  -  bv2)  , 


leaving 


1*1 


s2  +  bv2 


2 


for  the  remaining  term.  By  choosing  signal  locations  at  different  points  on 
the  characteristic,  so  that 

V  vi 

v2*  V2  ' 


we  can  make  this  term  vanish  by  a  suitable  choice  for  b.  Further  signal 
components  are  superfluous,  hence  the  multiple  target  hypothesis  always  allows 
us  to  make  Min  L2  ■  0. 

For  any  hypothesis  concerning  the  signal  structure,  we  can  define  the  log 
likelihood  function 

A  s  log  Max  f(*j,  s2)  -  -  j  Min  L2  -  log  (4*2)  , 
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where  Max  and  Min  refer  to  signal  paraaeters  contained  in  L2.  For  the 
noise- alone  hypothesis ,  HO,  there  are  no  signal  components,  hence  no 
paraaeters,  and  we  have 

A0  -  -  j  L„2  -log(4»2)  , 

where,  as  before, 

i*2  -  l*i|2  +  l*2l 2  • 

When  one  signal  component  is  allowed  (hypothesis  H^,  confined  to  the 
characteristic  trajectory  of  the  system,  we  have  found  that 

Min  L2  -  < |«! | 2  +  |*2|2)  sin2 (A/2)  , 

where  A  is  the  minimum  distance  attained  by  the  ML  estimate  of  location,  and 
hence  we  cab  wvlte 

Ai  ■  -  y  Lo2  sin2(A/2)  -  log(4w2)  . 

» 

Finally,  for  the  multiple-target  hypothesis,  H2,  we  have 
A2  -  -  log  (4w2) , 

since  we  have  shown  that  Min  L2  ■  0  in  this  case. 

In  pairwise  hypothesis  testing,  such  as  H1  verus  HQ ,  we  compare  the 
difference  of  the  log  likelihood  functions  to  a  threshold.  Thus,  we  accept  Hi 
over  HO  if 

A 1  “A®  >  v  . 


L-!i 


for  some  constraint  p,  which  will  be  chosen  to  meet  the  false  elan 
probability  requirements.  The  resulting  test ,  In  this  case,  has  already  been 
discussed.  To  generalise  the  testing  to  multiple  hypotheses,  we  Introduce  a 
constant,  p^,  for  each  hypothesis,  and  consider  the  numbers  -  p^,  for  k-0 , 
1  and  2.  If  the  largest  of  these  occur  for  k  -  1,  then  hypothesis  is 
chosen.  This  procedure  generalises  pairwise  testing  in  a  natural  way,  and 
(with  an  arbitrary  tie- breaking  rule)  leads  to  an  unambiguous  choice  of 
hypothesis  in  every  Instance.  This  procedure  can  also  be  viewed  as  a  special 
case  of  Bayes  hypothesis  testing  with  a  suitable  cost  matrix. 

Since  only  differences  of  the  p^  affect  the  decision  process,  we  can 
take  the  largest  of  the  numbers 

To  ■  “  y  Lo2  -  Po  » 

Tj  -  -  j  L02  sln2(A/2)  -  Pi  ,  and 

T2  “  *  M2  . 


% 


as  the  basis  for  decision.  Finally,  we  can  add  p£  to  each  of  these  quantities 
and  baae  decision  on  the  smallest  of  the  numbers 

to  5  Lo2  -  Co  , 


ti  s  Lq2  sln2(A/2)  -  Ci  ,  and 


tj  £  0. 


We  have  Introduced  new  threshold  constants,  CO-  2  (p2  -  pO)  and  Cl-2(p2-pl). 

A 

It  Is  interesting  that  L  2  and  A  are  sufficient  statistics  for  the 

°  2 

three-way  decision  process.  L.  is  a  measure  of  the  total  energy  contained 

A  V 

in  the  two  samples,  and  A  Is  a  measure  of  how  closely  the  data  ratio,  *^*2  > 
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resembles  an  Ideal  signal  ratio,  i.e.,  a  point  on  the  characteristic 
trajectory.  With  two  signals  present,  the  data  ratio  can  be  anything,  even  in 
the  absence  of  noise.  If  this  ratio  falls  on  the  characteristic  by  accident, 
the  data  will  look  like  a  single  target,  and  there  is  no  reason  to  expect  the 
receiver  to  decide  otherwise.  But  usually,  the  data  ratio  will  off  the 
characteristic  when  an  interferor  is  present,  and  A  will  allow  us  to  sense 
this,  even  though  the  received  energy  is  large.  In  a  simple  test  of  Hi  versus 
HO,  the  receiver  will  choose  Hi  if  the  received  energy  is  large  enough,  even 
the  estimated  location  is  far  from  the  system  trajectory. 

To  see  how  the  decision  mechanism  operates,  we  can  define  the  equivalent 
pair  of  statistics 

X  =  L0  cos (A/2)  , 

Y  =  L0  sin  ( A/2 )  ,  and 

plot  the  decision  regions  in  the  X-Y  plane.  Since  distances  between  points  on 
a  sphere  never  exceed  x,  both  X  and  Y  are  non-negative.  The  H0:H1  boundary  is 
given  by  to  -  tl  ,  or 

L02cos2(A/2)  -  C0  -  Ci  . 

Obviously,  CO  must  exceed  Cl  or  we  could  never  have  to  less  than  tl.  Hence 
the  required  boundary  is  the  line 


X  *  yj Cq  —  Cl  S  Xq  • 

Similarly,  the  H1:H2  boundary  is  the  line 


**  A"  !lo 


» 


and  Cl  la  positive  to  allow  tha  possibility  of  having  tl  <  t2.  Finally, 

the  H0:H2  boundary  Is  the  circular  arc  /x2  +  Y2  -  >/C0  -Vxo2  +  YO2  . 

The  three  boundaries  meet  at  the  point  (XO,  YO),  and  a  typical  situation 
would  be  as  shown  In  the  figure,  where  eadh  region  Is  labeled  by  the 
hypothesis  chosen  for  corresponding  data  points*  Only  those  segments  of  the 
boundaries  where  the  first  choice  changes  are  shown. 


In  general,  small  saaple  energy  reaults  In  the  "noise-alone"  decision, 
while  large  energy  Implies  either  H*  or  H2.  For  to  be  preferred,  the  angle 
A  must  be  small  (the  data  point  must  be  near  the  system  trajectory),  and  the 
larger  the  received  energy,  the  smaller  A  must  be.  The  parameters  XO  and  YO 
allow  some  control  over  the  performance  of  the  system,  as  characterised  by 
errors  such  as  false  alarms  and  Incorrect  editing. 

In  Section  7  we  deal  with  the  problem  of  Interference  In  another  way.  The 
Interfering  sources  are  treated  as  stationary  noise  sources,  and  it  Is  assumed 
that  the  resulting  noise  covariance  matrix  Is  known.  Maximum  likelihood 
estimation  is  then  performed  in  the  context  of  known,  non-white  noises.  That 
approach  makes  sense  when  one  has  an  opportunity  to  observe  the  noise  and  then 
attempt  the  estimation  of  target  location  before  the  noise  has  changed 
statistically.  The  treatment  in  the  present  section  Is  complimentary  to  this, 
dealing  with  intermittent  and  unpredictable  interference  In  a  simple  way,  by 
recognizing  and  discarding  contaminated  samples. 


3.  IDEAL  PHASE-COMPARISON  AND  IDEAL  AMPLITUDE-COMPARISON  MONOPULSE 

A  phase-comparison  aonopulse  system  is  essentially  an  interferometer, 
vhich  employs  two  identical  antennas  with  separated  phase  centers.  If  the 
complex  voltage  gain  of  either  antenna,  relative  to  its  own  phase  center,  is 
V,  then  the  corresponding  gains  relative  to  a  phase  center  midway  between  them 


can  be  written 


Vl(y)  -  V  e 


ikcosy 


V2(y)  -  V  e 


-ikcosy 


The  gain  function  V  depends  implicitly  on  the  actual  direction  of  arrival  of 
the  signal,  and  it  is  assumed  that  the  antenna  boresight  directions  are 
parallel.  Our  parameter  y  is  the  angle  between  the  signal  direction  of 
arrival  and  the  baseline  direction  established  by  the  line  between  the 
separate  antenna  phase  centers.  The  constant  k,  of  course,  is  ir  d/X  ,  where 
d  is  antenna  spacing  and  X  is  the  signal  carrier  wavelength. 

The  normalized  gain  functions  are 


v  (y) 
1 


1  ikcos  y  +  a 

m  P 


v  (y) 
2 


1  -ikcos  y  +  a 
6 


where  a  is  the  argument  of  V.  Only  y  can  be  measured  by  the  system,  of 
course,  and  the  gain  ratio  depends  on  this  directional  parameter  alone: 


vi(y) 

v2(y) 


_  2 ikcos  y 


The  characteristic  trajectory  of  this  system  is,  of  course,  the  equator  (or 
part  of  it)  on  the  unit  sphere: 

0  -  x/2 


2  k  cos  y 


If  d  exceeds  y,  Che  whole  equator  is  possible  and  y  Is  a  multivalued  function 
of  4.  As  mentioned  in  Section  2,  this  kind  of  ambiguity  cannot  be  resolved  by 
the  system,  and  Che  best  we  can  do  is  to  estimate  4*  By  restricting  our 
discussion  to  the  estimation  of  4,  we  are  dealing  with  all  phase-comparison 
sytem8,  and  the  final  translation  of  this  estimator  and  its  performnce  back  to 
the  parameter  y  will  not  be  carried  out.  He  also  choose  not  to  discuss  the 
case  d  <  y  explicitly,  and  define  the  "ideal  phase-comparison  system"  as  one 
whose  characteristic  trajectory  coincides  with  the  equator  of  the  unit 
sphere . 

For  the  ideal  phase-comparison  system,  we  obtain  immediately  the 
well-known  result  that  the  ML  estimator  of  4  is 

4  -  4  -  arg  (*x/*2)  -  «g  (*2**j>  * 

since  the  projection  of  a  point  on  the  sphere  to  its  equator  is  along  a 
meridian,  preserving  azimuthal  angle,  or  longitude.  Since  A  *  |n  -  •J  |  in 
this  case,  we  also  obtain  the  detection  criterion 

1  2  2 

2  (|zil  +  | Z2 |  )(1  +  »in  n)>  const. 

But 

2  2 

(|*il  +  U2I  )*in  n  ■  2 |  , 

and  therefore  detection  is  based  upon  the  equivalent  rule 

| zi |  +  | z 2 1  >  const. 

In  amplitude-comparison  monopulse,  two  antennas  with  different  patterns, 
but  sharing  a  common  phase  center,  are  employed.  In  standard  notation,  the 
gain  of  the  "sum  beam”  is  £,  and  that  of  the  "difference  beam”  is  A.  He  take 
Vi  ■  A  and  V2"  £  »  so  that 
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The  angle  y  Is  usually  azimuth,  and  Che  key  assumption  Is  that  Che  gain  ratio 
Is  a  real  function  of  y.  This  is  generally  true  in  practical 
amplitude-comparison  systems  within  the  sum  beam  main  lobe,  and  we  define  the 
"ideal  amplitude-comparion  system"  to  be  one  for  which  the  gain  ratio  assumes 
all  real  values  as  y  varies.  Ambiguities  are  again  possible,  but  the  only 
measurement  supported  by  the  system  is  an  estimation  of  the  value  of  this  real 
ratio.  Usually,  A  is  an  odd  function  of  y  and  £  is  even,  hence  we  assume  that 
y  ■  0  corresponds  to  a  zero  value  of  the  ratio. 

The  characteristic  trajectory  of  an  ideal  amplitude-comparison  system  is 
obviously  a  great  circle  through  the  poles  of  the  unit  sphere,  composed  of  the 
meridians  t  "  0  and  $  ■  w.  It  is  more  convenient  to  use  the  equivalent 
description  $  -  0,  allowing  0  to  range  from  -  ir  to  +  ir,  and  then 

ViCy)  -  tan  (8/2)  . 

V2(T) 

As  before,  we  study  amplitude-comparison  systems  in  general  by  dealing  only 
with  the  estimation  of  0.  The  properties  of  individual  antenna  systems  enter 
only  when  0  is  expressed  in  terms  of  y.  This  relation  is  usually  simple,  and 
near  boresight  is  often  well-approximated  by  a  linear  dependence  of  0  on  y. 

It  is  instructive  to  rederlve  the  known  results  for  amplitude-comparison 

A 

monopulse  from  our  formulation.  The  estimate  0  is,  of  course,  obtained  by 
dropping  a  perpendicular  from  the  data  point  (n  ,♦)  to  the  characteristic 
meridian.  By  elementary  spherical  trigonometry,  the  result  is  determined  from 
the  equation 

tan  §  -  tan  n  cos  ^ 


It  is  convenient  to  restrict  i|>  to  the  hemisphere 


*  ,  .  ^  * 

"  T  <  ♦  <  T  » 


and  let  n  range  over  [-ir,  w]  ,  in  which  case  it  will  be  found  that  0  is 
uniquely  determined  by  assigning  it  to  the  same  quadrant  as  n.  It  is  also 
seen  that  the  distance.  A,  to  the  characteristic  never  exceeds  ■§•  and  is 
determined  (using  the  Law  of  Sines)  by 

sin  A  -  |sin  n  sin  ^|. 

Amplitude-comparison  monopulse  is  usually  discussed  in  terms  of  the  gain 
ratio,  as  a  real  parameter  u,  and  the  data  ratio,  as  a  complex  variable  w. 
Then 


u  -  tan  (0/2) 
will  be  expressed  in  terms  of 

w  “  tan  (n/2)ei^  . 

This  connection  can  be  derived  by  noting  that 
tan  8  »  ■  - - 

1  -(5)2 

and  that 


2  Re(w) 

i  -H4 


2  tan  (n/2)cos  ♦ 
1  -tan*(n/2) 


tan  n  cos  t  . 


The  solutions  of  the  resulting  equation: 

“  ■  Re(w) 

1  -(5)2  I  -  jw|2 
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But 


Re2(w)  -  |  |w| 2  +  j  Re  (w2)  , 


and  hence  ( |w| 2—1) 2  +  4  Re2(w)  -  |w|4  +  1  +  2  Re  (w2)  -  |w2  +l|2  , 

which  gives 

G  -  M2  -1*  l»2  +  l\ 

±  1  Re(w) 


The  sign  ambiguity  Is  resolved  by  considering  the  case 
w  -  x  •  real. 

Then  u+  -  x  and  u-  -  -1/x  ,  and  only  the  positive  root  Is  consistent  with 

A 

the  constraint  of  keeping  8  and  n  in  the  sane  quadrant.  Thus 

G  -  t»2  +  *1  t  H2  -  2  Re  (w) 

2  Re  (w) 


|w^  +  l|  +  1  ■*  |  w|  * 


(1) 


the  desired  result 

As  for  the  detection  criterion,  we  note  that 


sinn 


1* 

e  • 


2  zlz*2 

1*1 1 2  +  I *2 1 2 


2  w 

TT-RT 


so  that 


sin  l  -  | slnn  sin*|  -  | 
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Then 


t  (1  +  |w|2)2  -4  Im2(w) 
cos  A  -  - 1 — 5 — 

1  +  |w2| 


1 1  4-  w2 1 
1  +  |w| 2  , 


as  a  consequence  of  the  identity 

Im2(w)  -  j  |w|2-  j  Re(w2)  , 

(2 ) 

and  we  obtain  again  the  known  result  for  the  detection  statistic: 

Y  (|zll2  +  lz2l2)d  +  cos  1 

*il2  +  I *al f 1  +  ' 

*  j- t |*i | z  +  U2I2  +  l*i2 

The  approximation  u  *  Re  (w)  is  often  used  in  applications  of 
amplitude-comparison  monopulse.  The  approximation  is  good  whenever  |w|«l 
(often  the  case  for  a  signal  near  boresight)  and-  it  is  exact,  as  we  have  seen, 
when  w  happens  to  be  real.  In  terras  of  the  sphere,  the  approximation  takes 
the  form 


+  z\\)  . 


tan(9/2)  »  tan  (n/2)  cos 41  , 

which  is  not  a  natural  one  to  make  in  this  context. 

The  analysis  given  in  this  section  is  aimed  at  developing  familiarity 
with  the  formulation  of  monopulse  problems  on  the  unit  sphere.  The 
idealizations  made  here,  of  two  special  great-circle  characteristic 
trajectories,  are  commonly  made,  in  some  form,  in  the  conventional  analyses  of 
these  problems.  They  are  often  acceptable  because  they  fail  significantly 
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only  outside  the  range  In  y  at  which  the  system  normally  operates.  One 
usually  does  not  expect  amplitude-comparison  monopulse  to  work  In  the  antenna 
backlobes,  for  example,  and  past  tracking  data  Is  almost  always  employed  In 
some  way  to  restrict  operation  to  a  well-behaved ,  main-lobe  portion  of  the 
coverage  which  also  permits  unambiguous  conversion  of  y  to  actual  signal 
direction. 


4.  TRANSFORMATION  OF  COORDINATES  IN  MONOPULSE  ANALYSIS 

It  has  been  known  for  a  long  time  that  the  analysis  of  an  amplitude- 
comparison  system  with  sum-beam  output  £  and  difference-beam  output  A  can  be 
transformed  into  the  analysis  of  a  phase-comparison  system  by  the  introduction 
of  the  new  quantities 


yr 


(£  ±  iA). 


The  signal  components  of  the  new  "coordinates"  are  equal  in  magnitude  and  the 
noise  components  are  again  independent  with  the  same  variances  as  before. 
Formulated  in  terms  of  S+  and  S_,  the  problem  is  mathematically  identical  to 
one  of  pure  phase-comparison  monopulse.  This  is  not,  of  course,  an  accident. 
In  this  Section  we  introduce  a  class  of  simple,  linear  transformations  of  the 
data  pair,  zj  and  z^  ,  and  investigate  the  induced  transformation  of 
coordinates  on  the  unit  sphere.  It  turns  out  that  data  transformations  which 
preserve  the  "white"  character  of  the  noise  samples  (independence  and 
equlvariance)  induce  rigid  rotations  of  the  sphere  and,  in  particular,  the 
transformation  from  (E,A)  to  (S+,  S_)  can  be  derived  as  the  one  necessary  to 
rotate  the  characteristic  trajectory  of  ideal  amplitude-comparison  monopulse 
into  the  equatorial  characteristic  of  its  pure-phase  counterpart. 


We  consider  the  samples,  z\  and  23,  to  be  components  of  a  vector,  z,  in 
a  two-dimensional  complex  space,  writing 


z 

-  1 

L  *2  J 


Similarly,  the  normalized  gains  and  the  noise  samples  are  components  of 
vectors  v  and  n,  and  we  have 


•  bv  +  n 


The  conjugate  row-vector  to  z  is  z  ,  given  by 


and  the  covariance  matrix  of  the  noise  vector  is  denoted 
M  ■  Enn*. 

We  have  modeled  this  noise  as  white,  each  complex  component  having  variance  2, 
so  that 

M  -  2  I2  , 

where  lg  denotes  the  2x2  unit  matrix. 

f  • 

Suppose  we  define  new  variables,  and  ,  as  linear  transforms  of  Zj 
and  Z£  by  writing 


where  T  is  a  non-singular  2x2  matrix  with  complex  components.  In  terms  of 
»  » 

z  ,  the  system  has  gain  vector  v  “Tv,  and  a  noise  covariance  matrix 

•  t  t 

M  -  TMT  -  2TT  , 

where  T*  is  the  complex  transpose  of  T.  Obviously,  if  we  choose  T  to  be 
unitary,  l.e.,  TT*  -  ,  than  the  new  noise  covariance  matrix  is  the  same  as 

the  old,  and  moreover  the  new  gain  vector  is  automatically  normalized. 

To  any  vector,  such  as  z,  there  corresponds  a  point  (?!,♦)  on  the  unit 
sphere,  determined  by  the  equation 

tan(n/2)e  ^  ■ 

z2 
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This  correspondence  is  many- to- one ,  since  z  can  be  multiplied  by  a  complex 
scalar  without  changing  n  and  The  transformed  vector,  *'■  Tz,  determines 
another  point  (n',  on  the  sphere,  through 


tan(n/2)e1*  '  - 

*2 


and  hence  the  transformation  T  induces  a  mapping  of  the  sphere  onto  Itself. 
This  mapping  is  given  explicitly  by  the  equation 


« 

tanCn^Je1* 


Tn  tanCn^Je1*  +  T^ 
Tjj  tanCn^Je1*  +  T22 


» 


in  terms  of  the  components  of  T.  This  relation  will  be  used  repeatedly  in 
this  study. 

Now  it  is  an  established  fact  that  when  T  is  unitary,  the  Induced 
transformation  on  the  sphere  is  a  simple  rotation.  This  equivalence  is  easily 
proved,  as  follows.  Let 

*1 

a  ’  [y 

by  a  vector  in  our  sample  space,  and  consider  the  inner  product 


If 


(a,z) 


*2*  *2 


tan(0/2)e1^ 


» 
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1*1 1  ♦  1*2! 


sin  (6/2) 


1*2 1 

1*1 l^+  l*2l 1 


2ai&2 


-  cos  (6/2)  ,  and 


•  slnBe  , 


exactly  like  the  relatione  used  in  Section  2.  Then,  writing  ||a||  ■  (a, a) 
2 

and  | | z | |  ■  (z,z)  for  the  noras  of  these  vectors,  we  coapute 


2  2  Z  2  2  *  A 

I  (a**)  |  “  |*i  |  |*i  |  +  |*2l  4*21  +2  **{*1  *2*1*2  } 

N 


Ia||2||*||2 


cose  1  -  cosn  ^  1  +  COS0  1  +  COST) 

5 - 2 —  + - 2 - 2 — 

+  Y  sln6  slnn  Refe1^  ^}J  . 

In  teras  of  A  ,  the  distance  on  the  surface  of  the  sphere  between  the  points 
(e,$)  and  (n,t)  ,  we  have 


iu,.)i2-  ii.|i2ii,na 


since 


cosA  ■  cose  cosn  +  sine  slnn  cos  ($  -  H»)  . 

Now  suppose  we  Introduce  a  unitary  transforaatlon,  V,  and  define  the  new 

1  1 

vectors  a*  “  Va  and  s'  “  Vs.  To  these  vectors  correspond  points  (e,  $ ,)  and 

•  1 

(n*  ♦)  on  the  sphere,  and  obviously 


I 


lu\.')l2.  Il.•ll^ll*'ll2  1  * f°“’ 

y 

where  A'  is  the  distance  between  (6',  +')  and  (n’»  ♦’).  But  a  unitary 
transforaation  leaves  inner  products  and  noras  unchanged,  and  we  can 
iaaediately  conclude  that  A'  ■  A  ,  so  the  transforaation  of  the  sphere  has 
preserved  the  distance  between  the  original  two  points.  Since  these  points 
were  arbitrary,  the  transforaation  aust  be  a  rotation. 

The  fact  of  this  equivalence  is  all  that  is  needed  in  this  entire 
analysis,  although  it  is  useful  to  know  exactly  which  unitary  aatrlx 
corresponds  to  a  prescribed  rotation.  A  slaple  derivation  of  this  detailed 
correspondence  appears  in  Appendix  A,  where  the  amplitude-comparison  /phase- 
coaparlson  transforaation  is  given  as  an  exaaple. 

From  an  analytical  point  of  view,  ideal  phase-coaparison  and  aaplitude- 
coaparlson  syateas  are  versions  of  the  same  basic  problea,  expressed  in 
different  coordinate  systeas.  Since  any  rotation  of  the  sphere  can  be 
effected  by  a  slaple  unitary  change  of  basis  in  the  saaple  data  space,  we  can 
say  that  all  aonopulse  systeas  are  equivalent  whose  characteristic 
trajectories  are  great  circles  on  the  sphere.  We  refer  to  these  as  "great 
circle  systeas",  and  say  we  are  working  in  "ideal  phase-coaparison  coordinates" 
when  we  have  rotated  the  characteristic  trajectory  onto  the  equator.  The 
important  portion  of  a  non-ideal  trajectory  of  a  real  system  can  be 
approximated  by  a  great  circle  arc,  and  then  rotated  onto  a  portion  of  the 
equator.  The  ideal  phase-comparison  estimator,  4  ■  4»  can  then  be  used  as 
an  approximation,  and  expressed  back  in  terms  of  the  original  sample  variables 
of  the  system. 

A  broader  class,  whose  Importance  will  appear  later.  Includes  all  systems 
whose  characteristic  trajectories  are  small  circles  on  the  sphere.  By 
rotation,  these  are  all  equivalent  to  systems  whose  characteristics  are 
circles  of  constant  polar  angle  (parallels  of  latitude),  say  6  -  6o.  It  is 
obvious  that  the  ML  estimator  of  4  for  such  a  system  is  still  4*4,  and  that 
the  minimum  distance  attained  by  this  estimate  is  A  »  |r»  -  6o|.  From  this 
last  it  is  easily  shown  that  the  detection  criterion  is 


1*1 |  sin(0o/2)  +  1*21  co8(0o/2)  >  threshold 


In  Part  II  it  is  shorn  that  the  accuracy  of  the  ML  estimator  of  0  for  this 
system  degrades  as  the  characteristic  trajectory  recedes  from  the  equator. 

For  this  and  other  reasons  it  appears  that  the  great  circle  trajectory  is  the 
best  that  can  be  attained  with  a  two  beam  monopulse  system. 

So  far  we  have  emphasised  unitary  transformations  of  the  data  vector 
because  they  leave  the  noise  statistically  invariant  and  preserve  the 
normalizations  of  the  gain  vector,  while  moving  the  characteristic  trajectory 
(whatever  it  may  be)  around  on  the  sphere  in  a  simple  and  useful  way. 
Non-unltary  transformations  are  also  important,  and  they  are  required  in  the 
whitening  operation  used  in  the  analysis  of  monopulse  in  Interference,  in 
Section  7.  These  transformations  sap  the  sphere  onto  Itself  in  a  more 
complicated  way,  but  the  mappings  are  still  conformal  and  circles  on  the 
sphere  are  still  transformed  into  circles,  although  without  the  preservation 
of  radius.  A  great  circle  trajectory  is  transformed  by  a  non-unltary,  noise- 
whitening  transformation  into  a  small  circle,  and  this  explains  the  appearance 
of  "small  circle  systems"  in  the  study.  The  basic  theory  of  non-unitary 
transformations,  as  applied  to  our  unit  sphere,  is  the  theory  of  linear 
fractional  transformations  of  the  complex  plane.  If  we  write  (*i/*2)  “  w  and 
(*l'/*2')  ■  w* ,  then  clearly  *’  -  T*  implies 


w' 


w  + 
w~+ 


t 


as  we  have  already  seen.  The  properties  of  these  transformations  are 
(3) 

well-known  ,  and  carry  over  onto  the  sphere  by  means  of  the  conformality 
and  circle-preserving  properties  of  the  stereographic  projection.  Some  of  the 
analysis  in  this  study,  particularly  in  Section  7,  is  motivated  by  the  known 
properties  of  linear  fractional  transformations.  The  exposition  here  is 
self-contained,  however,  with  proofs  of  statements  given  in  terms  of  spherical 
coordinates,  when  needed. 


5.  EXTENSIONS  TO  MULTIPLE-PULSE  ANGLE  ESTIMATION 

In  this  Section  we  make  two  extensions  o£  the  analysis,  each 
characterized  by  the  use  of  multiple  sample-pairs  as  inputs  to  the  estimation 
problem.  Our  first  example  is  rather  simple,  and  represents  a  class  of 
problems  exhibiting  signal  coherence  over  the  set  of  samples.  The  point  of 
the  example  is  to  illustrate  the  problem  of  the  simultaneous  estimation  of 
angle  and  other  signal  parameters,  such  as  Doppler  frequency.  Without 
additional  parameters,  nothing  new  would  be  Introduced  by  multiple  coherent 
samples,  which  is  a  situation  already  included  in  the  formulation  given,  since 
the  original  signal  waveform  was  quite  general. 

We  formulate  the  problem  directly  in  terms  of  the  sequence  of  sample 
pairs,  zi(n)  and  z2(n),  where  the  subscripts  refer  to  beams,  or  channels,  as 
before,  and  n  is  a  sample  index  running  from  1  to  N.  We  postulate  that  the 
signal  components  of  these  samples  have  the  form 


Ezk(n)  -  b  S(n)vk(y)  ;  (k  -  1,2)  , 


where  b  is  a  complex  amplitude  parameter  and  V}(y)  and  v2(y)  are  normalized 
beam  gains,  just  as  before.  The  new  feature  is  a  complex  signal  sequence, 
S(n),  presumed  to  depend  upon  one  or  more  implicit  parameters,  a,  and 
normalized  as  follows: 

* 

N  2 

I  |S(n)  |  -  1. 

n«I 


The  motivating  example  is  Doppler  modulation,  in  which  a  sequence  of  returned 
radar  pulses  is  passed  through  a  matched  filter,  as  described  in  Section  2, 
with  the  result  that  the  complex  output  amplitudes  are  modulated  by  a  Doppler 
sequence,  such  as 


S(n) 


n  o 
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In  all  cases  the  noise  is  assumed  to  be  white,  and  independent  and 
statistically  identical  from  sample  to  sample. 


Maximum  likelihood  estimation  of  all  the  unknown  parameters  immediately 
reduces  to  the  minimization  of 


-i.{ 


E  ijzxCn)  -  b  S(n)v1(y)|  +  |z2(n)  -  b  S(n)v2(y)|  }. 


|z}(n) |  +  |z2(n) |  +  | b  | 


I 


-  2  Re  b 


r(  n  *  *  *  *  ) 
<  I  Iv^yjS  (n)z1(n)  +  v2(y)S  (n)z2(n)]  % 


We  have  made  use  of  our  normalization  conventions,  and  we  now  Introduce  the 
definitions 


llzj^s  E  |zk(n)|2  , 
n-1 


ll*lll  +  1 1 *2 II 


» 


and 


u  * 

ffc  =  (S,zk)  =  E  S  (n)zk(n) 
n-1 

The  unknown  parameters  represented  by  a  are  now  implicit  in  f}and  f2. 
Substituting,  we  obtain 

I*  ■  L2  +  |b|2-  2  Re{b*lvx*(y)f i  +  v2*(y)  f2]}  , 

which  yields 

*  * 

bj  -  vx  fj  +  v2  f2 
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and 


2  2  ,  *  *  .2 
Min  L  -  L  -  |vx  £i  +  v2  f 2 I 
b  ® 

A 

The  ultimate  estimates  of  a  and  y  Are  to  be  substituted  into  bj  to  yield  b. 

2 

Note  that  our  expression  for  Min  L  is  exactly  the  same  as  the  one  derived 

b 

in  Section  2,  except  that  fj  and  f2  replace  and  *2.  Thus  we  could 
Interpret  the  zfc(n)  as  pairs  of  time  samples  of  the  incoming  signals  at  the 
antenna  terminals.  The  independent  noise  samples  would  then  be  direct  samples 
of  wide  band  white  noise.  The  f^  are,  of  course,  correlations  of  the  samples 
with  an  expected  waveform,  with  or  without  parameters  like  Doppler.  Then,  by 
a  somewhat  heuristic  limiting  process ,  one  could  conclude  that  matched 
filtering  was  a  good  thing  to  do  before  procedlng  with  angle  estimation.  But 
this  was  the  processing  postulated  in  Section  2 ,  and  the  point  of  the  present 
discussion  is  that  matched  filtering  is  a  proper  part  of  the  overall  ML 
estimation  problem,  starting  with  the  incident  time  waveforms. 

The  characteristic  trajectory  of  the  system  is  the  same  as  before, 
introduced  by  way  of  the  definition 

vi(y)  i* 

^  -  fa  (9/2)  .  . 

But  now  the  data  is  represented  by  the  point  (n»  ♦) ,  determined  by  the 
definition: 

fl(«) 

f^oT  "  tan  <n/2)e  v  , 

which  depends  upon  a,  yet  to  be  estimated.  If  A  the  arc  distance  between  the 
points  (n,f)  end  (6,4)  Are  the  unit  sphere,  then  we  have 

2  2  i  2  2 

Min  L  -  Lo-j<|fl|  +  |f2|  )(1  +  cos  A)  , 

since  the  expressions  are  the  same  as  in  Section  2  with  the  f^  replacing  the 
*k  • 
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With  an  arbitrary  system  trajectory,  the  estimates  of  y  and  a  will  be 
coupled  in  a  complicated  way,  since  A  will  be  the  minimum  distance  between  the 
characteristic  trajectory  and  the  curve  or  region  represented  by  the  points 
(n,  ♦).  However,  for  a  simple  trajectory,  the  estimates  can  be  separated  and 
performed  in  a  direct  sequential  way.  Using  the  small-circle  system  as  an 
example,  the  characteristic  trajectory  is  fixed  by  0  -  Q^and  for  given  a,  $  is 
estimated  as 

*  * 

♦  ■  ♦  ■  arg  <f2  fi)  , 

A 

and  A  *  |h  “  8  |  .  Finally,  a  is  estimated  by  maximizing  the  expression 
j  <|fi|Z+  |f2l  2)0  +  cos  A) 

-  5-  ( |f il 2  +  |f2|2)  +  J  cos  eo(|f2|2-  |fi|2>  +  sineo|flf2*| 

-  [sin  (0o/2)|f1(a)|  +  cos(eo/2)|f2(a)|]2  . 

For  an  ideal  phase-comparison  system,  this,  of  course,  reduces  to 
( I f 1 <a) |  +  |f2(o)|)2  . 

A 

The  final  estimate  of  4  is  the  value  of  $  evaluated  at  the  ML  estimate  of  a. 

In  our  second  example  there  are  no  other  parameters,  but  the  sample  pairs 
are  modeled  as  incoherent,  with  as  many  unknown  complex  amplitudes  as  there 
are  sample  pairs.  Specifically,  the  signal  components  are  taken  to  be 

E  ■  *k<n)  -  b(n)  vfc(Y)  , 

with  all  the  b(n)  unknown,  and  the  same  noise  model  as  in  the  previous 
example .  Then 


RBSH5S3 


>  ii  in  mi  w  n.m  fP  in'  «  »v  r**  ■*  ,  •  .»  v.  v 

.*4 .  *.*  h  .•••»■•  %  •  •>*.  >/  v“  •_*  *..  !*•  %  'An  . 


1 


i 


*  * . 


L  -  L  +  Z  |b(n) |  -2  Re  t  b  (n)  [v!  (Y)zi(n)  +  v2  (y)  *2(q)1 
°  n-1  n-1 


| 


and  the  amplitude  estimates  are 


bj(n)  -  v^CyJziCn)  +  v£  (y)  z2(n) 


The  estimator,  y  ,  will  minimize 


$ 


I 

> 

* 

s 

> 

*s 

4 


2  2  n  *  *  2 
Min  L  -  L  -  E  |vj  (y)zi(n)  +  v2  (Y)z2(n)|  ,  where  the 

0  n-1 

Min  is  taken  over  all  values  of  the  amplitude  parameters. 

He  expand  and  write 

Min  L2  -  -  |v!(y)|2|  z l  |2-|v2(y)f  I  z2|2  -2  Re  |v1*(y)v2(y)(z2,  Zj^) 


where 


(z2,z1)  =  Z  zl(n)z2  (n)  . 

n-1 

The  reduction  of  the  data  to  an  equivalent  point  on  the  sphere  Is  less 
straightforward  now,  since 

2  2  2 
|(z2,z1)|  #  1  z1  I  I  z2  I 

2  2  * 

except  as  a  very  special  case.  But  substituting  for  (v^j,  |v2|  and  vjv2 
In  terms  of  8  and  $  ,  according  to  the  formulas  introduced  In  Section  2,  we 
find 


Min 


L  -  L  -  A  I  z ,1  (1  -  CO80 )  -  i  I  z,l  (1  +  CO80) 
0  2  £  *• 

-  sin  0  Re  |e-i^(z2,zj)| 
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.a.V.I.  ,1.  rTm  iTi 


where 


j-  LQ2  -  j  (I  z2  I  “  »  zi  I  )  cose  -  8in8|(z2»2l>|co8(4>  -  ♦)  , 

4*  =  arg  (z2,  z0  • 

He  define  Q  and  n  by  the  equations 

2  2 

Q  COST)  s  I  Zj  I  -  I  *1  I 

Q  slnn  5  2 | (z2  »*i)  |  > 

and  then 

2  1  2  1 

Min  L  ■  4r  L  -  tt  Q  cosA  , 
tot 


where 


cosA  •  cosd  cosn  +  slnd  slnn  cos  ($  -  ♦) 


as  before.  Of  course, 

Vf  2  22  : 

(I  *2  ■  “  «  *i  ■  )  +  4 1 C*2 ,*i) | 


and 


tan(n/2)  e1* 


sin  n  e^ 
1  +  cosn 


2(*2 »  *1> 

~ - 2~  2 

Q  +  lz2l  -lljl 


2 

When  N-l,  Q  reduces  to  Lq  and  the  right  side  of  the  formula  above  becomes 
•Imply  Z}/z2,  as  It  mist. 

For  a  general  trajectory,  the  estimation  of  y  will  be  a  very  complicated 
function  of  the  data,  but  for  a  great-circle  or  small-circle  system,  using 

A 

phase-comparison  coordinates,  8  will  be  constant  and  $  given  by  the  relatively 
simple  expression 

A 

♦  -  ♦  -  arg  (z2 »*i)  . 


The  corresponding  detection  statistic  Is 


|  Q  cos(n  -  0Q) 

1  2  2  i  2 

*  j  (I  zj  I  +  I  *2  I  )  "*■  y  CO®  0Q(I  z2  *  “  *  *iO 
+  sin0o|(22»zl)| 

2  2  2  2 
-  sin  (0  /2 ) I  z,  I  +  cos  (0  /2)l  z2  I 
o  1  o 

+  2  sin  (0o/2)cos(0o/2)|(z2,*i)| 


The  ML  estimator  ij>  can  be  expressed  in  terns  of  the  Individual  estimators  for 
each  sample  pair,  namely 

-  l  *  \  (  *l<n> 

*n  =  arg  ^z2  (n)Zl(n)^  -  arg  j 

since  clearly 

♦  -  arg  |  I  |z1(n)z2(n)  |  e1*11 

(n*l 

A 

This  expression  shows  f  as  a  nice  mixture  of  coherent  and  incoherent 
processing.  The  individual  $  -  estimates  are  combined  in  an  essentially 
coherent  way,  as  a  weighted  sum,  which  is  appropriate  since  the  signal  ratio 
is  constant  for  all  the  pulses.  The  samples  appear  as  weights  in  an 
Incoherent  form,  since  signal  phase  is  unknown  and  variable  from  pulse  to 


6.  GENERALIZATION  OF  THE  NOISE  MODEL 

The  remaining  sections  of  Part  I  are  concerned  with  non-white  noise. 

The  natural  cause  of  correlated  noise  samples  will  be  external  sources  of 
deliberate  or  accidental  Interference.  Various  noise  models  are  discussed  in 
this  section,  and  it  is  shown  that  the  general  case  is  equivalent  to  a  mixture 
of  white  internal  noise  and  a  single  external  noise  source. 

For  the  general  noise  covariance  matrix. 

Mi  l  M12 

M2 i  M22 

the  diagonal  elements  oust  be  real  and  M21  must  equal  M12*,  since  M  is 
Hermltlan.  Moreover,  M  is  positive  definite,  which  (in  the  two-dimensional 
case)  is  equivalent  to  the  requirement  of  positive  trace  and  positive 
determinant.  If  we  write 


b  ■  d  cos8 
c  ■  d  sin© 
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and 


a  -  d  +  £  , 

with  d  and  £  non-negative.  Thus  one  general  fora  for  N  is 
£  +  d  (1  -  cos0)  d  sin0  e1* 

d  sine  e_i*  f  +  d  (1  +  cose) 

with  positive  f  and  d  and  unrestricted  values  of  6  and  4. 

Now  consider  an  external  source  of  noise,  at  azimuth  Yj  »  so  that  in  the 
absence  of  signal  the  total  effective  aodulation  at  the  antenna  terminals  is 
given  by 

\  -  Wj(t)  v^Yj)  +  Wk(t)  ;  k  -  1,  2  . 

Here  the  W^(t)  are  white  noise  Modulations  as  in  Section  2,  and  WJ(t)  is  an 
independent  white  noise  process,  with  spectral  density  N ..  : 

K  WjCtJW^t')  -  2  Nj  d(t-t'). 

Note  that  Wj  is  used  with  noraallsed  beam  gains. 

Following  Section  2 ,  we  assume  a  matched  f liter  system,  with  output 
samples 

*k  "  K  +  \ 

The  white  noise  components  are  unchanged ,  with 

Bl“ki2  “  2  , 

and  the  interference  components  are 

jk  -  G  vk(Yj)  /  s\t)Wj(t)dt  . 


It  follows  that  these  components  are  zero-mean ,  circular  Gaussian  samples  with 
covariance  matrix 


I 


e  jk  3» 


vk<^ J>vi  *YJ*  G  E  I  /  s 


2N 


r  V',j)vt*(YJ  > 

o 


as  a  result  of  our  filter  gain  normalization  convention.  The  location  of  the 
interferor  is  characterized  on  the  unit  sphere  by  the  point  (0j  $j)  ,  where 


Vl(YJ)  =  tan(0J/2)e  1<<>J 


V2(Yt> 


so  that 


ivl<Yj)|  -Id-  cos  0j)  , 


| v2CTj> 1 2  -  j  (1  +  cos  0 j)  ,  and 
Vi(Yj)v2(Yj)  -  j  sin  0 


In  terms  of  the  parameter 


j  =  i  i 

2  N  9 


the  covariance  matrix  of  the  total  noise  is 


M  -  2 


1  +  J  ( l-cosOj) 


J  sin  0j  e 


J  sin  0je 


1  +  J(1  +  CO80 j) , 


Comparison  with  the  general  form  shows  that  any  noise  model  can  be 
expressed  as  a  superposition  of  white  noise  and  an  external  source,  provided 
that  source  can  be  placed  anywhere  on  the  unit  sphere,  and  that  a  suitable 
scale  factor,  common  to  both  channels,  is  introduced.  This  scale  change  will 
only  affect  the  significance  of  the  signal  amplitude  parameter. 

It  is  not  illogical  to  place  the  interferor  off  the  characteristic 
trajectory  of  the  system,  since  that  trajectory  msy  well  be  an  idealization 
only  valid  in  the  system  main  lobe,  and  used  only  in  the  estimation  of  signal 
azimuth.  In  any  case,  as  we  shall  see  below,  the  resultant  of  a  number  of 
external  sources,  all  on  the  trajectory,  is  itself  off  the  trajectory  except 
in  a  special  case. 

The  significance  of  the  parameter  J,  essentially  a  ratio  of  intensity  of 
external  to  Internal  noise,  may  be  seen  by  noting  that 

e(|Ji|2+  lJ2|2) 

j  - - 2 - 2 - 

E  ( (n^ |  +  (na(  ) 


When  signal  components  are  added,  they  have  the  same  form  as  before,  namely 


where 


E  \  * b  • 


9 


and  Es  is  still  the  total  energy  collected  by  the  two  beams.  The  effective 
slgnal-to-noise  ratio  of  the  sample  pair,  when  both  signal  and  interferor  are 
present,  may  be  taken  to  be 


SNR 


|b«i|2+  lEg2f 
H1 1  +  M22 

E, 

2No+  NJ 

By  means  of  these  relations,  the  theory  can  be  applied  to  any  model  of  signal 
and  correlated  noise,  but  we  retain  the  notation  of  the  noise  covariance  last 
given,  in  terms  of  J,  0j  and  +j,  because  it  is  intuitively  useful.  Note  that 

Tr  M  -  4(1  +  J)  , 

and  that 

Det(M)  -  4(1  +  2J). 

It  is  instructive  to  examine  the  case  of  multiple  external  sources,  which 

are  characterized  by  locations  (0  ,  *  )  and  relative  intensities  J  ,  where  s 

8  8  8 

runs  from  1  to  S.  The  sources  are  assumed  to  be  independent,  and  the  total 
covariance  matrix,  Including  normalized  Internal  noise,  will  then  have 
components 


ibt 


4(1  +  J) 


Mu  -  2  |  1  +  E^  J8(l  -  cos  0g)  J 


M22  “  2 


(1  +  cos  0g) 


M12  -  2  E  J  sin  0  e  *■**  ,  and 

s-1  8  8 


m21  "  m12 


We  equate  this  to  the  matrix 
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1  +  J(l-C08  8j) 
J  sin  0je  ^ 


J  sin  Qj  e  i*J 
!  +  J(1  +  cos  8j) 


which  corresponds  to  our  standard  form  with  a  scale  factor  k  on  each  channel. 
The  multiple  sources  will  thus  correspond  to  an  equivalent  single  source, 
together  with  a  modified  Internal  white  noise  component.  The  connection 
equations  are 


k2J  sin  0j  e1**1  - 


E  J  sin  0  se1*8 
s-l  8  8 


2 

k  J  cos  8j 


S 

E  J  cos  0S 
s-l 


2  S 

k  Cl  +  J  )  -  1  +  E  J 

s-l  8 


and 


Equating  arguments  on  each  side  of  the  first  equation  yields  4Jf  and  the 

equality  of  the  remaining  magnitudes,  together  with  the  second  equation,  fixes 

2  2 

6j  and  the  product  k  J.  Finally,  k  is  found  from  the  last  equation. 

It  is  interesting  to  find  the  location  of  the  equivalent  external  source 
on  the  unit  spere.  Let  uj  be  a  unit  vector  whose  rectangular  components  are 
given  by 

<uj>*  -  cos  0j 

"  8in  0J  C08 

(uj)y  -  9j  0j  , 


and  let  u  ,  s-l,  ...8,  be  unit  vectors  corresponding  to  the  angles  (0S,+S). 
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We  equate  real  and  imaginary  parts  of  the  two  sides  of  our  first  equation 
above,  and  combine  these  two  with  the  next  equation,  supplying  appropriate 
unit  vectors,  to  find 


2 

k  J  u. 
— vJ 


S 

Z  J  u 
8“1 


Thus,  for  two  sources,  the  resultant  lies  on  the  great  circle  arc  joining  the 
respective  points  on  the  surface  of  the  sphere.  This  resultant  is  combined 
with  a  third  source,  and  so  on.  The  resultant  of  all  sources  lies  within  the 
smallest  convex  region  which  includes  the  individual  locations. 

With  a  group  of  external  sources,  each  located  on  the  system  trajectory, 
the  resultant  will  generally  be  off  that  trajectory.  The  one  exception  is 
the  great  circle  characteristic. 


7.  MAXIMUM  LIKELIHOOD  ESTIMATION  OF  ANGLE  IN  THE  PRESENCE  OF  INTERFERENCE 

We  hove  shown  in  Section  6  that  the  general  case  of  correlated  noise  in 
the  two  monopulse  channels  is  equivalent  to  a  nod el  having  white  noise 
components,  together  with  a  single  external  source  of  appropriate  strength 
and  location.  Specifically,  the  complex  samples  upon  which  estimation  will  be 
based  are  taken  to  be 


\  •  *k  +  4  +  “k  !  k  '  *• 2  > 


where 


8k  -  b  vk(Y)  , 


Jk  "  Vk  (V 


Here  b  is  an  unknown  complex  signal  amplitude,  y  represents  the  signal  angle 
(azimuth  or  elevation)  to  be  estimated,  yj  is  the  known  angle  of  the 
Interference,  and  j0  is  a  complex,  circular  Gaussian  variable  representing  the 
Interference  amplitude  after  filtering.  The  strength  of  this  Interference  has 
been  characterized  by  the  parameter  J,  in  such  a  way  that 


2 

•UJ  ■  -r  ■ 


As  usual,  the  white  noise  terms  nj  and  n2  are  independent,  complex  circular 
Gaussian  variables,  with 


2  2 

E|n}|  -  E|n2|  *  2. 


The  particular  definition  used  for  J  is  influenced  by  the  desire  to  keep  the 
formulas  which  appear  later  in  this  section  as  simple  as  possible.  Its 
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relation  to  signal- to-noise  ratio  la  discussed  in  Section  6.  In  this  section 
b  Is  an  unknown  nuisance  parameter ,  but  its  connection  with  signal  energy  and 
Internal  noise  level  given  by 


» 


is  important  in  the  analysis  of  estimator  performance.  With  this  model,  the 
total  noise  covariance  matrix  is 


M  -  2 


1  +  J(1  -cos0j) 
J  sinOje"1^ 


i*J 

J  sinflje 

1  +  J(1  +  cose.) 

J  — 


where 


/■< 


I 


i 


i*T  V1<YJ> 

t*n(9j/2)e 


In  order  to  solve  the  problem  of  estimating  signal  location  in  the 
presence  of  such  correlated  noise,  assumed  to  be  known,  we  seek  a 
transformation,  W,  with  the  property  that  in  the  transformed  coordinates, 

z'  «  Wz  , 

the  noise  is  white.  The  new  covariance  matrix  will  be 


M’ 


Ez 


(z’)1  -  W  M  Wf  , 
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* m  * k.  •  «  •«- 


and  we  require  Chat 


M’  -  2  I2  , 

so  that  estimation  problem  in  the  new  coordinates  will  be  of  exactly  the  type 
we  have  already  discussed.  The  transformation  matrix,  W,  cannot  be  unitary, 
since  a  unitary  transformation  would  leave  n  statistically  unchanged  while 
making  j  look  like  a  sample  vector  for  an  external  source  at  another  location. 
In  fact,  the  transformed  covariance  matrix  would  look  just  like  M,  with  the 
same  value  of  J  and  a  new  location  (8'j  ,  $'j  ),  fixed  by  the  effect  of  the 
corresponding  rotation  on  (8j  ,  ). 

Since  W  is  not  unitary,  we  do  not  know,  in  general,  what  the 
characteristic  trajectory  of  the  system  will  be  in  the  new,  whitened, 
coordinates.  Moreover,  there  are  many  choices  for  W,  indicating  the 
possibility  of  many  different  characteristics  in  the  new  coordinate  system. 
Fortunately,  all  the  whitened  coordinate  systems  are  simply  related  and,  in 
fact,  one  can  be  transformed  Into  any  other  by  means  of  a  unitary 
transformation.  Thus,  the  various  "whitened"  characteristics  on  the  sphere 
differ  by  simple  rotations.  To  see  this,  note  that  if  W  is  a  whitening 
transformation,  then 

2  w"1(Wt)"1  -  M  , 

and  hence 

V^W  -  2M_1  . 

If  V  is  another  whitening  transformation,  then 

V*V  -  WfW  , 
from  which  it  follows  that 


which  implies  that 

(V  W  -i)-1  -  (V  tf"V  . 

Thus  V  W  *  is  unitary  and  since 

z"  S  V  z  -  (V  W_1)  s'  , 

the  result  is  proved.  This  fact  allows  us  to  choose  W  in  such  a  way  that 
the  characteristic  trajectory  is  changed  in  a  controlled  and  desirable  way. 
When  the  estimator  is  then  transformed  back  to  the  original  coordinates  and 
expressed  in  terms  of  the  original  data,  the  result  will  be  Independent  of  the 
choice  made  for  W,  so  long  as  W  whitens  the  noise. 

Perhaps  the  simplest  example  of  non-white  noise  is  described  by  a 
diagonal  covariance  matrix.  The  noise  components  are  independent,  but  have 
unequal  variance.  In  terms  of  our  general  covariance  matrix,  the  equivalent 
external  source  is  at  one  of  the  poles  of  the  sphere.  Such  a  source  would  be 
in  a  null  of  one  of  the  antenna  beams.  The  simplest  whitening  matrix  for  this 
case  is  also  diagonal. 


with  the  scaling  factors,  a  and  b,  chosen  to  adjust  each  noise  component  to 
the  standard  variance,  namely  2. 

The  mapping  of  points  (6,  $)  into  new  points  (6*,  $')  on  the  sphere  is 
very  simple  in  this  case.  The  general  transformation  equation  reduces  to 


tan  (eWe1*'  -  |  tan(e/2)e^  , 

for  this  situation.  The  longitudes  of  points  are  unchanged  ($•■»  $) ,  while  the 
latitudes  are  distorted  by  the  non-linear  transformation  of  the  polar  angles: 


tan(8'/2)  ■  tan  (0/2)  . 

D 


If,  for  example,  a  >  b,  then  all  points  are  moved  toward  the  south  pole,  and 
points  with  small  values  of  0  move  more  than  points  farther  from  the  north 
pole. 

A  meridian,  such  as  the  characteristic  trajectory  of  an  ideal 
amplitude-comparison  system,  is  invariant  as  a  curve,  but  individual  points 
are  transformed  according  to  the  above  equation.  A  parallel  of  latitude  is 
displaced  to  a  new  latitude,  but  the  relationship  of  points  before  and  after 
the  transformation  is  particularly  simple,  since  longitude  is  preserved. 
Actually,  this  transformation  (indeed,  any  linear  transformation  of  the 
coordinates)  takes  circles  on  the  sphere  into  other  circles,  but  we  can 
arrange  things  so  that  these  general  properties  are  not  required  to  solve  the 
estimation  problem. 

To  deal  with  the  general  coveriance  matrix,  we  proceed  indirectly,  asking 
first  what  transformations  (beside  rotations)  leave  the  equator  of  the  sphere 
invariant  as  a  curve.  We  Intend  to  deal  only  with  the  case  of  an  Ideal  phase- 
comparison  system  in  general  noise,  since  any  system  with  a  circular 
characteristic  can  be  reduced  to  this  case  by  a  rotation  followed  by  a 
diagonal  transformation,  and  the  resulting  noise  covariance  matrix  can  then  be 
represented  in  our  general  form.  The  effect  of  a  general  transformation, 


on  an  equatorial  point,  (ir/2  ,$),  Is  given  by  the  equation 


tan(8  72)e 


a  e*^  +  b 


c  e1^  +  d 


If  we  choose  c  ■  b*  and  d  ■  a*  ,  the  right  side  can  be  written 


a  e 


* 

b  e 


+ 


+ 


9 


which  obviously  has  magnitude  unity.  Thus 


is  an  equator-preserving  transformation,  and  it  is  not  hard  to  show  that  the 
most  general  transformation  with  this  property  has  the  form  of  W,  followed  by 
a  rotation  about  the  z-axis  of  the  sphere.  This  last  rotation  is  of  no  help 
in  whitening  the  noise,  hence  we  see  what  can  be  done  with  the  form  for  W 
given  above.  Since  a  matrix  to  whiten  M  satisfies  W*W  -  2M”*  ,  we  find  out 
what  W  can  do  by  computing 

1*1 2  +  M2  2  a*b  1 


2  ab  *  |a|2  +  |b|2 
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1  +  J  (1  +  cosQj)  -J  sinQJe 


WfW  -  2M~ 


1  +  2  J 


-J  sin0je"i0J  1  +  J(1  -  cosOj) 


Let  us  define 


and  put 


-  s  -  tan  (0/2 )e*^ 
a 


1*1 2  +  |b|2  -  1  , 


since  a  constant  factor  can  bt.  absorbed  in  the  quantities  d+  and  d-  .  Then 

2  2  _  2  2  _ 

| a |  -  cos  (0/2),  | b |  -  sin  (0/2)  , 


*  -  i* 

2  a  b  ■  -  sin  0e  *  ,  and 


d+  +  d_  +  (d+  ~  d_)cos  0 


-(d+  +  d-)sln  0e*0 


W*W  -  - 
2 


“(d+  +  d_)sin  0e  d+  +  d_  -  (d+  -  d.)cos  0 


Comparison  with  the  desired  form  shows  that  we  must  choose 


♦  ■  ♦j  » 
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and  (equating  the  traces). 


1  +  J 


7  (d+  +  d_)  - - 


1  +  2J 


Finally,  we  must  have 


sin  6 


J  sin  6j 
1  +  J 


7  (d+  -  d_) 


1  +  2J 


J  COS  0J 


cos  6 


We  find  8  from  the  equation  for  sin  0,  taking  the  solution  in  the  first 
quadrant,  and  solve  for  the  diagonal  factors: 


Notice  that 


2  - 
cos  8 


1  +  J  (  J  cos  0j  ) 

1  +  2J  ^  (1  +  J)  cos  0  l 

(J  sin  0j  \2  1  +  2J  +  J 

1  +  J  /  (1  +  J 


2  2  a 

cos  8j_ 


so  that 


1  +  2J  +  J  cos  8j 


V”  * 


The  diagonal  elements  are  equal  In  thi8  matrix,  and  the  same  will  be  true  of 
its  inverse,  hence  W  can  only  whiten  noise  covariance  matrices  which 
correspond  to  external  sources  on  the  equator  (to  make  cos0j  >0).  For  a 
perfect  phase-comparison  system  this  would  be  the  case,  but  as  discussed  in 
Section  6,  we  wish  to  analyze  the  model  having  an  ideal  phase-comparison 
trajectory  but  arbitrary  noise. 

Although  a  matrix  like  W  cannot  whiten  the  general  noise  covariance,  it 
turns  out  that  such  a  W  can  diagonalize  the  general  M.  To  put  it  another  way 
M  can  be  whitened  by  a  sequence  of  two  transformations,  the  first  of  which 
preserves  the  equator,  while  the  second  is  an  appropriate  diagonal  matrix.  To 
see  this,  we  write  the  general  such  two-stage  transformation  matrix  in  the 
form 


where  d+  and  d_  are  additional  parameters,  which  will  be  real  and  positive 
(complex  factors  can  be  reduced  to  this  case  by  z-axls  rotations).  Again  we 
compute  W^W,  and  find 

d+|a|  2  +  d_  |b| 2  (d++d_)a*b 

W*W  - 

.(d++d_)ab*  d+  |b| 2  +  d_  |a| 2 

Now  we  have  enough  generality  to  match  any  noise  covariance  matrix,  and  we 
require  that 


(J  COB  Qj  I 

- -  I  <  1 

(1  +  J)  cos  0  • 

so  that  d+  and  d-  are  positive* 

The  whitening  properties  of  V  depend  only  on  the  relative  phase  of  a  and 
b,  and  we  choose  to  make  a  real,  completing  the  specification  of  parameters  in 
W.  Altogether,  we  have 


The  first  stage  of  this  transformation  preserves  the  equator;  the  second  moves 
it  to  the  circle  parallel  to  the  equator,  a  distance  0’  from  the  north  pole, 
where 


yw 


/♦•  -  ♦A  i  + 

V7  /  "  77 


tan  (6/2) 
tan  (6/2) 


tan  I  — 


P) 


The  point  4  “  ♦j  on  the  equator,  and  Indeed  the  entire  meridian  containing  the 
source  location.  Is  Invariant  to  the  whitening  transformation,  thanks  to  our 
choice  of  real  a.  The  factor  connecting  the  tangents  can  be  simplified  by 
squaring  It,  with  the  result  that 


1  +  tan  (1/2) 

1  -  tan  (5/2) 


1  +  sin  6 


1  -  sin  6 


In  terms  of  the  original  noise  parameters,  the  mapping  of  points  (t/2,+)  on 
the  original  characteristic  to  their  images,  (O',*'),  on  the  characteristic  In 
whitened  coordinates,  Is  described  by  the  equations 


cosfl'  -  - 


J  cos  6j 


1  +  2  J  +  J2  cos2  0j 


/♦'  -  \  / 1  +  J  (1  + 

\  2  /  V  l  +  j  (l  - 


sin  9j)  /♦  -  4>j 

-  tan  I  -  -  - 

sin  0j)  \  z 


)  • 


A  signal  vector,  bv,  is  changed  by  this  transformation  to  b  W  v,  which  we 


write  as 


b'  v*  ■  b  W  v  . 

The  change  in  the  amplitude  parameter  le  needed  since  W  is  not  unitary.  The 
component  ratio  vj/v2  defines  a  point  on  the  original  trajectory: 


e1*  , 


while  v^'  /v2'  defines  the  image  point,  as  discussed  above.  The  new  vector, 
v' ,  is  normalised,  hence  the  new  signal  amplitude  satisfies 

|b'|2  -  |b|2  j  |Wnvi  +  W12v2|2  +  IW21VJ  +  W22  v2  |2j  . 


We  substitute  for  the  components  of  W  and  note  that 


|v2|2  -  1/2  ,2  V!  v2*  «  e  , 


with  the  results 


,  i  /  j  co.8j\ r  i  i 

l«n*i  +  -  tttt  +  -z 


sine  cos(f-fj) 


2  /  1  J  cos6j  \  f  1  l 

|W2l  V|+ W22  V2|  -^TTTT  1  +  J  “  Tose  y  [  2  "  2 


■in  e  CO. (4  -  fj) 


Altogether,  we  find 


1 

[  1 

I)1/2 

b*  -  b 

1  1  +  2J 

1  +  J  -  J  sin  0j  cos  (4  -  ♦j)|> 

• 

which  shows  that  the  new  signal  amplitude  depends  on  the  original  signal 
location.  This  adds  no  coapllcation  to  the  estimation  problem,  since  b*  Is 
Just  a  nuisance  parameter,  but  the  relation  of  b'  to  b  Is  required  to  evaluate 
the  performance  of  the  ML  estimator  In  this  case,  which  Is  done  In  Part  II. 

Now  that  we  have  a  suitable  whitening  transformation,  we  can  complete  the 
analysis  by  finding  the  ML  estimator  of  signal  location  In  the  new  coordinates 
and  transforming  It  back  to  the  original  set.  For  a  single  pulse,  the  data 
vector,  z.  Is  transformed  to  z’  -  Wz,  and  the  ML  estimator  on  the  whitened 
sphere  Is 


♦  '  -  arg  UV*'  l  *  » 

since  the  characteristic  Is  a  parallel  of  latitude.  In  terms  of  the  original 
data  components,  we  have 

erg  +  *22*2)*  0*n*l  +  W12*2>} 

-  arg  ||  -  sln(0/2)e  i*J  zj  +  cos(9/2)z2j  |cos(0/2)z1  -  •in(0/2)e1^Jz2j j  , 


since  positive  real  factors  do  not  affect  the  argument  of  the  quantity  on  the 
right.  Continuing,  we  get 


«•* 

* 


k> 


£ 
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♦  '  -  arg  |  *1*2  cos  (5/2)  +  *1  *2  sin  (8/2 


*  2  - 


^  2  2  *♦  J  | 

-  sin(8/2)cos(8/2)  | z  ^ |  +  |z2|  «  f 

arg  ^e1^*1  |  (1  +  cos  8)ziz2*e  +  (1  -  cos8)zi*z2  e1^ 

-  sin  8  (  |  z  2 1  +  | z2 1  ) J  | 


Finally,  we  factor  out  the  positive  quantity,  sine,  to  obtain 


ifj 


"  /  _  *  "tj  .  ft 

♦*  -  4j  ■  arg  J  cot(8/2)z1s2  e  +  tan  (8/2)z1  z2  e 

9 

-( 1*1 12  +  t *2 1 2 >  | 


On  the  original  sphere,  the  estiaator,  $  ,  is  given  by  the  Inverse  of  the 
whitening  transforaation: 


tan 


-  sin8 


tan 


sln8 


(^) 


The  source  longitude,  tj,  enters  as  a  sort  of  reference  angle  in  these 
expressions,  which  is  a  consequence  of  the  fact  that  the  source  aeridian  is 
Invariant  to  the  whitening  transforaation.  It  is  interesting  that  the  other 
noise  paraaeters,  J  and  8j,  enter  these  equations  only  through  8  ,  so  that  the 
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A 


final  estimator,  $,  for  the  general  case  Is  identical  in  form  to  the  estimator 
for  a  special  case  having  the  source  on  the  equator  at  the  same  longitude  and 
having  an  equivalent  intensity,  to  maintain  the  value  of  0.  The  performance 

A 

of  $  will,  of  course,  depend  expllclty  on  J  and  9j. 

A 

In  the  single-pulse  case  the  formula  for  can  be  written  so  that  it 
depends  on  the  data  only  through  the  ratio  zj/z2.  However,  in  its  present 
form,  this  expression  is  easily  generalized  to  the  case  of  multiple  pulses 
with  Independent  amplitudes,  as  in  the  second  problem  considered  in  Section  5. 
He  assume  the  noise  is  stationary  over  the  sequence  of  pulses,  so  the 

A 

whitening  transformation  is  the  same  for  every  pulse.  Then  Is  given  by 


♦  '  -  arg 


(  °  *<  ) 
j  I  z2  (n)  zj*  (n)  ) 


Repeating  all  the  steps  of  the  derivation  just  given,  we  easily  find  that 


♦  -  -  arg  |  cot (6/2) e  Z  z2*(n)zi(n)  +  tan(6/2)ei^J  I  z^(n)z2(n) 

I  n*l  n-1 


-  Z  (|*l(n)|  +  |z2(n) |  ) 

n-1  1 


In  the  notation  of  Section  5, 


I  z2  (n)z1(n)  -  (z2,  Zj) 

n-1 


*k(n) |  - 


•V2 


k  -  1  ;  2 


which  allows  us  to  write 

♦'  -  ♦j  -  arg  |  cot  (6/2)(*2.zi)  «  _1*J 

+  tan  (0/2)  (*2,*i)*e  -  I s x I  -  1*2*  }  • 

A  A 

The  transformation  from  4  to  the  desired  4  is,  of  course,  the  same  as  in  the 
single-pulse  problem. 

These  relations  provide  a  formal  solution  to  the  multiple-pulse 
estimation  problem  in  correlated  noise.  However,  a  more  succinct  result  is 
possible  by  following  a  slightly  different  line  of  reasoning.  Instead  of 
locating  the  estimate  as  a  point  on  the  system  trajectory  in  the  whitened 
coordinates  and  then  transforming  it  back  to  the  original  sphere  (according  to 
the  equations  just  derived),  it  proves  useful  to  consider  the  meridian  on  the 

A 

whitened  sphere  which  passes  through  the  characteristic  at  the  longitude,  4  , 
of  the  estimate.  This  meridian  maps  back  onto  the  original  sphere  as  a 
circle,  which  must  pass  through  the  equator  (the  original  characteristic)  at 

A 

the  longitude  4>  .  By  pursuing  the  course,  we  obtain  a  simple,  explicit 
* 

formula  for  t  ,  and  also  a  geometrical  picture  of  the  estimation  process  as  a 
projection  on  the  original  sphere. 

Without  loss  of  generality  we  can  take  ■  0,  since  the  original  sphere 
can  always  be  rotated  to  make  this  true,  and  the  meridian  corresponding  to 
(and  ir  +  $j)  was  preserved  by  the  whitening  transformation.  Next,  we 
introduce  the  notation 

2(*2»  *l) 

r  s  — — — — 

2  2 

Isxl  +  1*2*  , 

A 

and  notice  that  the  equation  for  +'  can  now  be  written 

♦  '  -  arg  {cot(0/2)  r  +  tan  (0/2)r*  -  2}  . 
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Thus  r  is  a  sufficient  statistic,  in  the  original  coordinates,  for  the 
estimation  of  signal  location.  The  meridian  through  the  ML  estimator  on  the 
whitened  sphere  is  the  set  of  points 

w'  ■  tan(0'/2)e^ 

for  all  values  of  8*.  Replacing  tan  (8'/2)  by  a  real  parameter  k,  and 
Introducing  the  definition 

y  =  tan  (8/2)  , 

we  can  write 


w'  -  k  -  r  +  yr  -2 


As  k  takes  all  positive  arei  negative  real  values,  w'  traces  out  the  desired 
meridian.  The  image  of  w'  on  the  original  sphere  is  w,  determined  by  the 
inverse  of  the  transformation  W.  It  is  easy  to  show  that 


Det  (W)  -  (1  +  2J) 


-  1/2 


and  that,  in  the  general  case, 


cos(8/2) 


sln(8/2)e 


/ Jcos8j 

1  +  J  -  - —  0 

CO80 


sin(0/2)e 


cos(0/2) 


0  \/l  +  J  + 


JcosOj 


We  substitute  in  the  Inverse  transformation 
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Wll1  W  +  vlz'1 


*21  *'  +  W22 


putting  ■  0  and  using  the  notation  d+  and  d_  Introduced  before.  The  result 


/ ST  cos(0/2)w'  +  /3+  sin0/2) 

w  " - 

/3l  Sln(0/2)w'  +  cos(0/2) 

Dividing  numerator  and  denominator  by  /d~cos(072)  and  recalling  our 
definition  for  y,  we  find 


1/2 

(d_/d+)  '  w'  +  y 
(d_/d+)1/2  m'  +  1 


Now  we  substitute  for  w'  ,  and  absorb  the  quantity  (d_/d+)  ’  into  the 
free  parameter  k  (without  changing  notation),  to  find 


k  (-  r  +  yr  -  2  +  y 
Uk  /'j  r  +  yr*  -  l\  +  1 


This  Is  a  parametric  representation  of  the  curve  on  the  original  sphere 
which  corresponds  to  the  meridian  of  the  estimator  on  the  whitened  sphere.  For 
the  moment,  we  put 
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so  ChaC 


1  * 

-r  +  yr  -  2  =  G  , 


kG  +  y 

ykG  +  T  * 


and  solve  as  follows: 


yw  -  1 


,  “  -  w  • 

k  w  -  ii 

The  statement  that  the  right  side  of  this  equation  is  real  is  the  equation  of 
our  curve.  It  may  be  written 

*  *  * 

(w  -  u)  ( yw  -  1)G  -  (w  -  y)  (yw  -  1)  G  , 


y(G  -  G*)  (jwj  +  1)  +  (G*-  y  G)  w  -  (G  -  y  G*)  w  -  0. 


G  -  G*  «  (  -jj  -  y)  (r  -  r*)  -  2i(  ^  -  y)  Im(r)  , 


2  *  i  3  2 

G-yG  -  (—  -  y  )  r  -  2  (1-y) 


■  y  ( 


+  y)  r  -  2 


«v  .-.V, /.  .  .  •*.  -\ 


hence  our  equation  simplifies  to 


2i  Im(r)(|w|2  +  1)  +  U)r*  -  2j  w 

-[(*  +  “MW*  ■  °- 


Recalling  the  definition  of  p,  we  find  that 


-jj-  +  p  ■  cot(0/2)  +  tan(0/2) 


2 

sin0 


and  therefore 


|w|  +  1  +  i 


_  *  * 

(r  -  sin0)  w  -  (r  -  sln0)w 
sinl  Im(r) 


-  0  . 


It  is  easy  to  put  this  equation  in  the  form 


which  represents  a  circle  in  the  w  -  plane  and  also  on  the  sphere,  after 
substituting 


w  »  tan^^Je1^  . 

However,  we  can  better  describe  this  circle  by  identifying  certain  points 
which  must  be  on  it.  If  w  is  real,  our  equation  becomes 

2  2w 

w  +  1  -  - —  ■  0  , 

sini 
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which  is  satisfied  by  the  values 

w+  “  tan(6/2)  ,  and 

w_  “  cot(0/2)  . 

These  are  the  values  w+  ■  p  and  w_  -  1/p  ,  which  correspond  to  w"  ■  0  and 
w'  ■  ~  ,  the  poles  of  the  whitened  sphere.  The  image  points,  w±  ,  lie  at 
equal  distance  from  the  poles  of  the  original  sphere,  on  the  meridian  4  ■  0. 
This  Is  the  meridian  of  the  interferor,  and  w  ^.represents  the  point  (8,  <t>j) 
in  general,  while  w_  represents  (ir-0,  4>j)»  A  line  through  w+  and  w-  is 
parallel  to  the  north  -  south  axis  of  the  sphere,  and  planes  through  this  line 
cut  the  sphere  in  all  possible  circles  containing  w+  and  w_  .  Our  equation 
describes  one  of  these  circles  which  passes  through  the  equator  at  the  point 

A 

The  circle  also  passes  through  a  point  which  depends  only  on  the  data 
(note  that  w+  and  w_  depend  only  on  the  noise).  To  find  this  point,  we  must 
pick  w  proportional  to  r,  since  that  will  cause  sin0  to  drop  out  of  the 
equation.  Putting  w  -  kr,  with  k  real,  we  have 

2,  ,2 

k  |  r  |  +  1  -  2k  -  0  , 

or 


k 


1  ± 


9 


hence  the  points. 


lie  on  the  circle.  They  are  equidistant  from  the  equator,  at  the  same 

A 

longitude.  Thus  t  °&y  he  determined  as  the  equatorial  Intercept  of  the  circle 
through  Wf  and  w_  which  passes  through  R+  and  R_.  If  the  Interference 
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goes  away,  i.e.,  J  ♦  0,  then  9  ♦  0  and  the  points  w+  and  w_  approach  the 

“  A 

poles.  Circles  through  w+  and  w_  then  become  meridians,  and  4  Is  simply 
arg(r)  -  arg(s2»cl)>  the  standard  white  noise  solution.  With  Interference 
present  the  lines  of  projection  (the  circles  through  w+  and  w-  )  become 

A 

distorted,  moving  $  away  from  the  value  It  would  have  In  white  noise. 

A 

To  get  an  explicit  expression  for  4  »  we  have  only  to  put 


w  ■  e 


to  find  where  the  circle  Intercepts  the  equator.  Substitution  yields 


(r  -  sln6)e  ^  -  (r*-  sin8)e  ■  21  sin©  Im(r) 


a  =  arg(r  -  sin 9)  , 


sln(a  -  $)  ■  sln0 


Im(r) 


| r  -  sinl | 


But  Im(r)  -  Im(r  -  slnO),  and  hence 


sln(o  -  ♦)  •  sln0  slna  , 


♦  ■  a  -  sin  (sln0  slna) 


■  arg(r  -  sln0)  -  sin" 


^  (  slni  Im(r)  | 

\  |r*  7 
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As  noted  before,  $  ♦  arg(r)  as  sin0  ♦  0,  which  naans  that  the  sin  on  the 
right  side  of  the  above  equation  must  be  taken  in  the  range  (—■ ir/2  ,  ir/2).  In 
the  general  case,  with  interference  at  longitude  $j,  the  result  is 

a  _ 

sin(a-$)  ■  sin0  sin(a  -$j)  , 

with 

a  =  arg(r-  sin0e  . 

In  the  single-pulse  case,  we  have 

2z2  zl  2(*/j*2) 

r  -  -  -  - 

2  2  2 

1*1 1  +  I *2 1  1  +  1*1 /*2 1 

and  it  follows  that 

R+  -  z\/  Z2  ■  t»n  , 

R_  -  z2*/*i*  ■  cotCn/ZJe1^  . 

This  simple  result,  that  the  projecting  circle  pass  through  the  data  point 
(n,^)  itself  in  the  one-pulse  case,  corresponds  to  the  fact  that,  on  the 
whitened  sphere,  the  estimator  lies  on  the  meridian  of  the  data  points 
*1 ,/*2,» 

It  will  be  noted  that  our  solution  of  the  correlated-noise  problem  is 
based  on  what  might  be  considered  a  fortunate  guess  for  the  whitening 
transformation.  This  approach  was  used  because  it  minimizes  the  amount  of 
tedious  algebraic  computation  required.  Another  method,  which  involves  a 
direct  construction  of  the  whitening  transformation,  depends  upon  the  detailed 


correspondence  between  rotations  of  the  sphere  and  unitary  transf oraations , 
which  is  developed  in  Appendix  A.  Using  this  correspondence,  one  can  write 
the  whitening  transformation  as  a  sequence  of  three  transformations,  as 
follows.  The  first  is  a  rotation  which  moves  the  location  of  the  interferor 
to  the  north  pole.  The  noise  ,covar lance  is  then  diagonal  and  the 
characteristic  is  a  great  circle  whose  center  is  known.  The  second  transform 
is  a  simple  diagonal  whitening  operation  which  changes  the  characteristic  to  a 
small  circle  whose  center  and  radius  can  be  directly  evaluated.  The  third 
transformation  is  another  rotation,  which  moves  the  center  of  the 
characteristic  circle  to  a  pole,  leaving  us  with  standard  white  noise  and  a 
trajectory  which  is  parallel  of  latitude.  It  is  possible  (but  not  simple)  to 
show  that  the  result  is  exactly  the  same  as  the  transformation,  W,  used  above. 

Still  another  method  is  to  begin  with  a  triangular  transformation  matrix 
which  removes  from  *2  its  projection  on  tj,  and  equalizes  the  resulting  noise 
variances.  This  operation  leaves  the  noise  white  and  the  characteristic 
circular.  By  keeping  track  of  what  has  happened  to  this  circle,  one  can 
rotate  it  into  position  as  a  parallel  of  latitude.  The  result  is  again  W,  and 
in  all  cases  the  meridian  of  the  Interferor  remains  invariant. 
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8.  NON-DETERMINI STIC  SIGNAL  MODELS 


In  all  the  cases  so  far  considered,  the  signals  have  been  modeled 
deterministically,  but  with  unknown  parameters  according  to  various 
hypotheses.  Instead,  we  can  assume  that  the  signal  waveform  is  a  random 
process,  such  as  white  noise,  or  simply  postulate  that  the  signal  components 
of  the  sample  pair  (or  pairs)  are  random  variables  with  a  corresponding  moment 
matrix.  We  have  already  shwon  (in  Section  6)  that  the  moment  matrix  of  one 
sample  pair  takes  the  following  form,  when  the  only  input  is  an  external 
source  of  white  noise: 


M  -  2 


1  +  J(I  -  cosflj) 
J  sinOJ  e~10J 


J  sin6J  ei8j 
1  +  J(1  +  CO80j) 


The  angles  8j  and  locate  the  source  on  the  unit  sphere,  and  we  now  assume 
that  this  point  lies  on  the  characteristic  trajectory  of  the  system.  This 
source  location  is  to  be  estimated,  and  we  allow  the  system  trajectory  to  be 
completely  general.  The  parameter  J  is  the  source-to-internal  noise  spectral 
density  ratio: 


o 

and  its  significance  was  discussed  in  Section  6.  The  samples  have  been 
normalised  so  that  M  -  2  I2  in  the  absence  of  any  external  source,  but  the 
model  for  M  is  completely  Independent  of  the  particular  form  of  the  filters 
used  in  the  system.  In  another  application,  one  might  wish  to  estimate  the 
azimuth  of  an  external  source  with  a  signal  waveform  different  from  that  for 
which  the  receiver  channels  are  matched.  This  signal  could  be  in  the  same 
band  but  utilize  a  different  spreading  code,  for  example.  The  response  of  the 
matched  filters  to  a  random  code  from  this  source  could  be  modeled  so  that  the 
samples  are  approximately  Gaussian  random  variables,  and  a  moment  matrix  of 
the  form  M  would  result,  with  an  appropriate  value  of  J. 


Some  of  Che  oldest  work  on  monopulse  has  used  Che  random  model,  and  we 
show  here  that  the  ML  estimator  of  source  location  Is  Identical  to  that 
obtained  with  deterministic  models,  regardless  of  the  nature  of  the 
characteristic  trajectory  of  the  system. 

With  this  random  model,  the  pdf  of  the  pair  of  samples,  and  Z2,  is 


-(ztM^z) 


f(s)  "  -  e 

ir2Det(M) 


where 


“  [z*]  and  Zt  “  [Zl**  *2*]  . 


as  before.  Note  that  f(z)  is  the  joint  pdf  of  four  real  Gaussian  variables. 
We  already  know  that 

Det(M)  «  4(1  +  2J)  , 


and  that 


2(1  +  2J) 


1  +  J(1  +  cos0 j)  -  J  sinOj  e*^ 


-  J  sin0j  e1^  1  +  J(1  -  cos6j) 


Substitution  yields 


logf(z)  -  ~ 


__i_j 

2(1  +  2J)  ( 


11  +  J(l  +  cosSjJJlzil2  +  [1  +  J(1  -  CO80j) ] | z2 | 2 


-  2J  sinfij  Re  |ziz2  e  !<>g(l  +  2J)  -  log(4ir2). 
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A 


With  our  standard  notation  for  the  data  ratio: 


—  ■  tan  (n/2)*1^  , 

*2 

we  can  use  the  expressions  found  in  Section  2  for  |zi|2  ,  | *2 1 2  *1*2*  » 

with  the  result 


i  1*1 12  +  l*2hf  /  ,  1 

logf(z)  -  -  j - mHi — I1  +  J  “  J(cos0jcosn  +  slnejsinncosdj-  ♦))! 


-  log(l  +  2J)  -  log(4x2). 

Again,  only  the  distance,  A  ,  between  the  points  (8j  ,  tj)  and  (n ,♦)  on  the 
unit  sphere  enters  this  expression,  and 


logf(z)  -  -  £ 


j  )*ll2+l*2l: 


1  +  2  J 


(1  +  J  -  JCO8A) 


-  log(l  +  2J)  -  log(4*2)  . 

The  ML  estimator  of  source  location  Is  the  same  as  before,  namely  that 
point  on  the  system  trajectory  nearest  the  data  point,  (n,  ♦).  If  A  Is  the 
resulting  minimum  distance,  then 


1  l*ll2  +  1*2 12  , 

Max  log  f(z)  -  -  j  - py  2  j —  |l  +  2Jsin2 

8J,  l 


-  log(l  +  2 J)  -  log(4x2)  . 


The  estimate  of  J  is  then  found  by  straightforward  differentiation: 
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J  -  t  |T  (l*il2  +  l*2l2)co®2  (f)"  i 

or 

"j'no  |r<l«il2  +  l*2|2w  (f)- 1 1 

For  an  ideal  phase-comparison  system ,  we  have  0^  ■  w/2  ,  ,  and 

A  “  ir/2  -  n  .  Then 


and 


cos* 


({)-  —i 


cosA  m  1  +  slnn 


“  No|  T  +  l*2|2)<1  ♦  •I*1*)  "  *J 
/ 

■  »0 1  i  <l*i I2  +  M2>  +  j  1*1*2* I  - 

a  ' 

■  H0 1  i  < 1*1 1  +  1*2! )2  ■ 


We  recall  that  sj  and  *2  normalised  filter  outputs ,  using  a  convenient 

choice  of  filter  gain.  In  terms  of  the  quantities 


\s  I  S*(t)*k(t)dt  , 


which  are  the  actual  outputs  of  filters  matched  to  S(t),  we  have 


This  analysis  is  easily  generalised  to  the  case  of  N  sample  pairs, 
assumed  independent  of  each  other,  and  statistically  identical  to  the  single 
pair  already  modeled.  If  the  nth  sample  pair  is  denoted  zi(n)  and  Z2(n),  as 
before,  then  the  logarithm  of  the  joint  pdf  of  these  samples  is  simply 

log  f  -  -  2TF+TJT  |  11  +  J  (1  +  cose^Jlzjl2  +  [I  +  J(1  -  coaeJ)]iz2i2 

-  2  Jsln8j  Re  je  (z2,zi) 

-  Nlog  (1  +  2J)  -  Nlog  (4ir2)  . 

In  this  expression  we  are  using  the  vector  space  notation  Introduced  in 
Section  5.  As  in  that  section,  we  define  the  effective  data  point,  (n,  ♦), 
by  the  equations 

Q  sinn  e1^  *  2  (z2,z1) 

Q  cosn  ■  lz2l2  -  Izil2  , 

with 

!Z  2)1/2 

(lz2l2  -  lzll2)  +  4 | Cz2,z1) j  l 

In  terms  of  these  quantise,  we  can  write 

log  f  -  -  2(1  +  2J)  |  ^  +  J^[,zl|2  +  lz2,Zl  "  JcosA 

-  Nlog(l  +  2J)  -  Nlog(4x2)  , 

where  A  is  the  distance  from  (6j,  to  the  point  (n,  ♦)  on  the  unit 
sphere. 


Again,  the  ML  estimator  of  source  location  is  the  same  as  In  the 
deterministic  case  discussed  before,  where  every  sample  pair  had  an  unknown 
amplitude  and  phase  factor.  If  we  define 

P  =  Izil2  +  |*2|2  # 

we  can  express  the  interference-level  estimator  in  the  form 
*  N 

N.  -  —  (P  +  Q  cos  A)  -  N  , 

J  4N  ° 

A 

where  A  is  the  minimum  distance  to  the  system  trajectory  attained  by  the  ML 
estimate  of  source  location.  For  an  ideal  phase-comparison  system,  we  again 

A 

have  cosA  ■  sinn  ,  and  then 

Nj  “  Nq  jjL  |*1»2  +  |*2|2  +  2|(z2,z1)|  -i 

a  direct  generalization  of  the  single-pair  result. 

Since  we  have  modeled  the  external  source  waveform  as  random  noise,  we 

can  interpret  this  source  as  an  interferor  instead  of  a  desired  signal,  and 

# 

make  use «of  the  above  analysis  to  estimate  its  location.  However,  in  this 
application  we  may  wish  to  remove  the  constraint-  that  the  source  lies  on  the 
system  trajectory,  since  the  source  may  be  a  composite  of  many  separate 
lnterferors.  In  Section  6  we  showed  that  the  resultant  may  be  off  the 
characteristic  in  this  case,  and  also  that  the  multiple  sources  affect  the 
level  of  the  apparent  white  noise  components. 

With  these  considerations  in  mind  we  postulate  that  the  noise  covariance 
has  the  more  general  form 

1  +  J(1  -  COS0.) 

2  J 
M  -  2k 

J  sinfij  e"1^ 


J  sineje1**1 
1  +  J(1  +  cosGj) 
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where  k  is  an  additional  parameter  corresponding  to  an  unknown  white  noise 
component.  This  is  exactly  the  form  used  in  Section  6  to  represent  the  total 
noise  equivalent  to  a  collection  of  external  sources  and  white  Internal  noise. 
We  assume  that  the  data  consists  of  N  sample  pairs  and  estimate  (0^, 
first,  then  J,  and  finally  k.  The  Intermediate  results  can  be  used  if  the 
subsequent  parameters  are  presumed  known.  When  all  the  parameters  are 
unknown,  M  becomes  an  arbitrary  positive-definite  matrix,  and  our  final  result 
reproduces  the  known  ML  estimator  of  M,  namely  the  sample  covariance  matrix. 

The  addition  of  the  factor  k2  is  a  simple  matter,  and  the  logarithm  of 
the  joint  pdf  becomes 

log  f  -  -  - - - |  (1  +  J) [I*! I2  +  l*2l2]  +  J  cose.llxxl2  -  I *2 * 2 1 

2k2(l+2J)  I  J 


-  2  J  sinOj  Rele”1*-1  (*2,*i)j|“  N  log[k**(l  +  2J)]  -  N  log(4ir2). 


With  the  Introduction  of  the  data  parameters  (n ,  ip) ,  P  and  Q,  as  before,  we 
obtain 

log  f  -  -  - i -  )(1  +  J)  P  -  JQ  cosA 

2k2 (1  +  2J)  { 

-  N  logk** ( 1  +  2J)  -  N  log(4v2)  . 


The  source  location,  (0j,  4j) »  Is  now  unconstrained,  hence  we  make  the 
eatlmates 


which  corresponds  to  A  ■  0.  With  this  choice,  we  have 

/ 
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log  f 


[P  +  J  (P  -  Q) J 


2k2(l  +  2J) 


-N  logk4(l  +  2J)  -  N  log<4*2) 


- 2_±_2 - N  iog(l  +  2J)-  ^  -  N  logk*»  -  N  log(4»2) 

4k2(l  +  2J)  4k2 


Differentiating  with  respect  to  (1  +  2J),  we  find  that  the  estimator  of 
J  satisfies  the  equation 


1  +  2J 


P  +  Q 
4Nk2 


which  leaves 


log  f  -  -  N  - 


N  -  N  log  /--*  -  N  logk2  -  N  log(4v2). 

\  4N  J  4k2 


If  we  stop  at  this  point,  i.e.,  treat  k  as  known,  the  estimated  covariance 
matrix  Is  rather  complicated ,  with  elements 


M11  “ 


+k2+  (t± 

4N  Q  V  41 


^2  -  ULfi  +  k2  + 

4N 


l*2l2  “  »*1»2 


(H?-) 


*  * 


M12  -  M21 


-  2 


<*2 


-  .•) 
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rv  .■% 


If,  however,  we  continue,  k  is  estimated  by  maximizing  the  expression 


-  "  N  logk2  , 

4  k* 


with  the  result 

k2  -  Ljl -a 

4N 

A 

When  this  estimate  Is  substituted  in  M,  the  complicated  quantity  Q  disappears 
from  the  formulas  and ,  substituting  for  P ,  we  obtain  the  expected  result 


M  - 


1 

N 


izi*2  (*2  »*l) 

(Z},Z2)  IZ2I2 


APPENDIX  A 


THE  CORRESPONDENCE  BETWEEN  UNITARY  TRANSFORMATION 
IN  THE  SAMPLE  SPACE  AND  ROTATIONS  OF  THE  UNIT  SPHERE 

In  Section  4  it  was  shown  that  a  unitary  matrix,  V,  viewed  as  a  linear 
transformation  in  a  two-dimensional  complex  vector  space,  corresponds  to  a 
rotation  of  the  unit  sphere.  The  mapping  of  a  vector 

*1: 

*  "  .Z2. 

to  a  point  (0  ,$)  on  the  sphere  is  expressed  by  the  relation 

Z1  i* 

—  ■  tan  (0/2)  e1?  . 

*2 

Any  linear  transformation  matrix  can  be  multiplied  by  a  complex  scalar 
constant  without  affecting  the  transformation  induced  on  the  sphere,  hence 

4  ^ 

e  V  and  V  correspond  to  the  same  rotation.  In  this  case  the  scalar  factor 
has  magnitude  unity  to  preserve  the  unitary  character  of  the  transformation. 
The  correspondence  of  transformations  can  be  made  unique  by  requiring  that 
Det(V)  ■  1,  which  we  now  assume. 

To  obtain  the  specific  correspondence  in  question  we  consider,  first, 
infinitesimal  transformations.  We  write 

V  -  I2  +  i  e  H  , 


and  find  that 

W*  -  I2  +  i  e  (H  -  Hr)  , 

to  first  order  in  e.  Thus  H  must  be  Hermitlan,  and  also,  since 
Det  (V)  -  1  +  i  e  Tr  (H)  , 
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i*i  .*t.> 


r_.  -> 


co  first  order,  Che  crace  of  H  must  be  zero.  Since  e  is  included  as  an 
explicit  factor,  the  elements  of  H  can  be  normalized  in  some  way,  and  it 
proves  convenient  to  use  the  following  parametrlzation  for  V: 


where 


V  -  I2  +  (i/2)e  S(]j)  , 


S<h)  s 


V3  “(w l  +  iP2) 

-(pi  -  ip2)  -  U 3 


pj2  +  p22  +  p32  =  1  . 


The  numbers  ui.  p2  and  M3  are  real,  and  may  be  interpreted  as  the  components 
of  a  unit  vector  u  in  a  real,  three-dimensional  space.  Obviously, 


Tr(S)  -  0 

and  also,  we  find  that 


2 

S  - 


The  induced  transformation  on  the  sphere  takes  (0,4)  into  (6 ’,4’), 

where 


V j i tan  (0/2 Je1^  +  V12 

tan  (O’^Je1*  -  - rr - 

V21tan  (0/2 )e  ’*  +  V22 


We  write  6'  *>  0  +  A6,  4> *  m  4>  +  A$  and  expand  the  left  side  to  first-order 
quantities.  We  also  substitute  for  the  components  of  V  and  develop  the  right 
side  of  the  above  equation  to  first  order  in  e.  The  result  is 

tan(0/2)  +  sec2(0/2)A0  +  i  tan(0/2)A«j>  +  ... 

■  tan(0/2)  +  iep3tan(0/2)  -  (ie/2)(ui  +  iy2)e 

+  (ie/2)(yi  -  ip2)tan2(0/2)ei^  +  ... 

or  A0  ■  e  C U2c°84>  “  Pisin$) 

A $  “  e  [ IJ3  -  (y^cos^  +  y2sin<(i)cot  • 

Now  if  we  write 

i* 

y3  -  co80q  ,  yj  +  iy2  ■  sin0o  e  TO  , 

id-.itifying  the  coordinate  axes  of  the  three-space  in  which  ]i  is  defined  with 
those  of  the  space  in  which  the  unit  sphere  is  embedded ,  we  have 

A0  -  e  8in0Q  sin($Q  -  $)  ,  and 

sin0A$  ■  e  [cos0Q  sin0  -  sin0o  cos0  cos($Q  -  $)] 

These  relations  show  that  the  point 

0*0  |  4  *  A 

o  T  To 

is  unchanged  by  the  transformation,  and  that  its  antipode 


is  also  invariant.  Since  V  corresponds  to  a  rotation,  these  points  lie  on  Its 

axis,  which  is  to  say  that  y  Itself  defines  the  axis  of  rotation  of  the 

sphere.  A  plane,  normal  to  y  and  passing  through  the  center  of  the  sphere 

defines  a  great  circle,  which  is  left  invariant  by  the  rotation.  At  the 

longitude  $  ■  +  *  ,  this  circle  reaches  its  highest  latitude,  being 

tangent  there  to  the  circle  0  ■  —  -  0  .  At  this  vertex,  we  find 

2 

A6  -  0  ,  and 

sin0  A<f>  -  e  sin  (0  +  0  )  ,  or 

CO80  A<k  •  e  , 
o 

which  shows  that  the  angle  of  rotation  corresponding  to  V  is  e  itself. 

Now  let  V  (y  ,  x)  be  the  unitary  matrix  corresponding  to  a  rotation 
through  a  finite  angle,  x>  ln  the  positive  sense,  about  the  axis  y  .  Then 

V(y,  dX)  -  I2  +  (i/2)dX  S(y)  . 

The  sequence  of  transformations  described  by  the  product  V(y,  dx)  V(jj  ,x)  then 
corresponds  to  a  rotation  through  x  about  y,  followed  by  a  rotation  through 
dx  about  the  same  axis.  The  resulting  rotation,  through  x  +  dx  about  y 
corresponds  to  V(ji  ,  x  +  dx)  ,  hence 

V(jJ  ,  X  +  <*X)  "  V(]j  ,  dx)  V(H  ,x) 

-  V(u  ,X)  +  U/2)dX  S(u)  V(H,x)  . 

Thus  V  (y  ,x)  satisfies  the  differential  equation 

j  V(M  »X)  -  (i/2)  S(H)  V(H  ,X) 
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with  V(u  ,0)  -  Ij 


The  8olutioo  is  simply 


V (jj  .  X)  ■  .<1/2)*S<1!)  .  I2 

*  S  il  °[S("^ 

n-o 


But  S2n  ■  I2  a-id  S2n  +  *  ■  S,  hence 


V(jj,  x)  ■  cos(x/2)  I2  +  i  Bin  (x/2)  S(w)  , 
or,  in  matrix  form 


V<H,  X) 


cos(x/2)  +  iP3  sin(x/2) 

(-ipj-  u2)  *in(x/2) 


(-iyi  +  y2)  sin(x/2) 

% 

cos(x/2)  -  in 3  sin(x/2) 


He  can  now  use  this  correspondence  to  derive  the  transformation  of  the 
data  vector  which  will  convert  an  ideal  amplitude-comparison  system  into  an 
ideal  phase-comparison  one  in  the  new  coordinates.  First,  the  sphere  must  be 
rotated  through  90°  about  the  X  axis,  in  order  to  carry  the  amplitude- 
comparison  trajectory  (which  lies  in  the  XZ  plane)  into  the  equatorial  plane. 
In  the  process,  the  north  pole,  which  corresponded  to  the  data  point  (A/E)~0, 
is  moved  to  the  intersection  of  the  negative  Y  axis  with  the  sphere.  A  second 
rotation  through  90°,  this  time  about  the  new  Z  axis,  will  preserve  the 
trajectory  while  moving  this  reference  point  to  the  nominal  reference,  6~tr/2, 
<(>*0,  for  an  ideal  phase-comparison  system.  From  the  general  form,  a  rotation 
about  the  X  axis  corresponds  to 


Vi(x> 


”cos(x/2) 
_-i  sin(x/2) 


-i  sin(x/2) 
cos(x/2) 


and  a  Z-axis  rotation  Is  given  by 

,ix/2) 

V3(x) 


- o  — 0 

r  ei*/2>  o  *i 
"  L  o  .-‘x'd  . 


The  complete  transformation  we  require  is  V$(— )  V^(— )  ,  or 

2  2 


i»/4 

e 


V  - 


0 


The  new  coordinates  are  then 


Z'x-  e_lir/4  1  (Z2  +  iZj) 


Z'o"  e“i,r/4 


—  (Z2  -  iZx>  . 

h 


-lff/4 

The  phase  factor,  e  ,  is  completely  harmless  and  can  be  dropped,  which 
completes  the  derivation  since  l\  was  the  difference  channel,  Z2  the  sum 
channel  sample  in  our  amplitude-comparison  formulation. 
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